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Experimental studies of superconductivity properties of fullerides are briefly reviewed. Theoretical 
calculations of the electron-phonon coupling, in particular for the intramolecular phonons, are 
discussed extensively. The calculations are compared with coupling constants deduced from a 
number of different experimental techniques. It is discussed why the A3C6O are riot Mott-Hubbard 
insulators, in spite of the large Coulomb interaction. Estimates of the Coulomb pseudopotential 
/i* , describing the effect of the Coulomb repulsion on the superconductivity, as well as possible 
electronic mechanisms for the superconductivity are reviewed. The calculation of various properties 
within the Migdal-Eliashberg theory and attempts to go beyond this theory are described. 
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I. INTRODUCTION 

The fullerenes (Cgo, C70,...) have attracted much inter- 
est since their discovery by Kroto et al. (1985), not least 



because of their appealing, symmetric shape (see Fig. 1). 
This interest increased dramatically when Kratschmer 
et al. (1990) discovered how to produce Cgo in large 
enough quantities to be able to make solids of a size 
which allowed traditional solid state experiments. Very 
soon Haddon et al. (1991) found that intercalation of 
alkali metal atoms in solid Ceo leads to metallic behav- 
ior. Shortly afterwards it was found that some of these 
alkali-doped Cgo compounds are superconducting with a 
transition temperature T c which is only surpassed by the 
cuprates (Hebard et al., 1991; Rosseinsky et al., 1991; 
Holczer et al., 1991a; Tanigaki et al., 1991; Fleming et 
al, 1991). Thus T c is 33 K for RbCs 2 C 60 (Tanigaki et 
al, 1991) and for CS3C60 under pressure T c = 40 K has 
been reported (Palstra et al, 1995). The great interest in 
the superconductivity of the alkali-doped Cgo compounds 
is in particular due to these systems being a completely 
new class of superconductors, the large value of T c and 
the question whether or not such a large value of T c can 
be caused by the coupling to phonons alone. There has 
therefore been a great effort over the last five years to 
characterize and understand both the normal state and 
superconducting properties of fullerides. The present re- 
view deals only with Ceo compounds, since no other su- 
perconducting fullerides are known at present. 




FIG. 1. The C 60 molecule. 



1 




FIG. 2. Some of the subbands around the Fermi energy for solid Ceo in the Fm3 structure. The bands at about -0.5 
eV are the h u bands, which are occupied in solid Ceo, and the bands at about 1.5 eV are the t\ u bands, which become 
populated in A n C6o- The bands around -1.5 eV results from the overlapping h g and g g bands. (From Erwin (1993).) 



Ceo is the most symmetric molecule in the sense that 
its point group (icosahedral) with 120 symmetry opera- 
tions is the largest point group of the known molecules. 
It has the shape of a soccer ball. The 60 carbon atoms are 
all equivalent and form 12 pentagons and 20 hexagons. 

The Ceo molecules condense to a solid of weakly bound 
molecules. While the shortest separation between two 
atoms on the same molecule is about 1.4 A, the shortest 
separation between two atoms on different molecules is 
about 3.1 A. The fullerites are therefore molecular solids, 
where many of the molecular properties essentially sur- 
vive in the solid. The discrete levels of a free Ceo molecule 
are only weakly broadened in the solid, leading to a set of 
essentially nonoverlapping bands with a width of about 
1/2 eV as is illustrated in Fig. 2. The system therefore 
has two very different energy scales, the intramolecular 
{Ei ~ 30 eV) and the intermolecular (W ~ 1/2 eV) en- 
ergy scales. For an undoped solid, the h u band is 
full and the t lu band is empty, and this system is there- 
fore a band insulator. When solid Cgo is doped by alkali 
atoms, the alkali atoms donate about one electron each 
to the t\ u band (Erwin and Pederson, 1991; Saito and 
Oshiyama, 1991b; Martins and Troullier, 1992; Satpa- 
thy et al, 1992). Since the t\ u band can take six elec- 
trons, it is half-full for A 3 Ceo (A=K, Rb), which is a 
metal (Haddon et al, 1991). A4C60 is, however, an in- 
sulator (Murphy et al, 1992; Kiefl et al., 1992; Bcnning 
et al, 1993) although the t\ u band is only partly filled 



and should be a metal according to band theory (Erwin, 
1993; Erwin and Bruder, 1994). A 6 C6o is a band insula- 
tor. Recently, ACgo has attracted much interest because 
of its rich and interesting phase diagram and because of 
its possibly quasi-one dimensional character (Stephens 
et al, 1994; Chauvet et al, 1994, Erwin et al, 1995). 
TDAE-Ceo has a ferromagnetic-like transition at T = 16 
K (Allemand et al, 1991). If the fulleridcs are considered 
as organic compounds, they have the highest supercon- 
ducting and (ferro) magnetic transition temperatures for 
organic compounds. 

In A3C60 (A=K, Rb) the Ceo molecules are located on 
a fee lattice. At low temperatures, the Ceo molecules take 
one of two likely orientations in an essentially random 
way (Stephens et al, 1991). This oricntational disorder 
has a substantial effect on some of the electronic prop- 
erties (Gclfand and Lu, 1992ab; 1993; Mele and Erwin, 
1994), but has been neglected in most calculations. A 
certain short-range correlation of the orientations is ex- 
pected on theoretical grounds (Gunnarsson et al, 1991; 
Mazin et al, 1993a, Erwin and Mele, 1994) and has been 
observed experimentally (Teslic et al, 1995). The alkali 
atoms are located in the tetrahedral and octahedral holes. 
In Rb NMR of Rb 3 Ceo one would then expect to observe 
two lines. At low temperatures, however, three lines are 
seen (Walstcdt et al, 1993; Zimmer et al, 1993), raising 
some questions about our understanding of the precise 
structure. 
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The effective Coulomb interaction between two elec- 
trons on a Ceo molecule in a solid is about 1-1^ eV (Lof 
et al, 1992; Briihwiler et al, 1992; Pederson and Quong, 
1992; Antropov et al, 1992). The long-ranged Coulomb 
interaction in the A3C60 compounds leads to a charge 
carrier plasmon, due to the oscillations of the t\ u elec- 
trons. This plasmon has the energy lu p i ~ 1/2 eV, an 
intermediate coupling constant (g/aj p i) 2 <~ 1 (Antropov 
et al, 1993b; Knupfcr et al, 1993) and unusal dispersion 
(Gunnarsson et al, 1996) and broadening (Liechtenstein 
et al, 1996). 

In systems where the superconductivity is driven by 
the electron-phonon interaction, the transition tempera- 
ture T c can be calculated in the Eliashberg (1960) theory. 
A simple estimate of T c can be obtained from the McMil- 
lan (1968) formula 



L2 



exp[ 



1.04(1 + A) 
A-/i*(l + 0.62A) J 



(1) 



where ui n is a typical phonon frequency (logarithmic av- 
erage), A is the electron-phonon coupling, and n* is the 
Coulomb pseudopotential, which describes the effects of 
the repulsive Coulomb interaction. To estimate T c we 
then in particular need to consider A and fi* . 

Cqo has intramolecular vibrations (phonons) with en- 
ergies up to ujph ^0.2 eV. Only phonons with A g or H g 
symmetry couple to the t\ u electrons. It was very early 
proposed that these phonons drive the superconductiv- 
ity (Varma et al, 1991; Schlutcr et al, 1992ab; Mazin et 
al., 1992). Estimates of the electron-phonon interaction 
A ~ I — 1 fall in the right range to explain the experimen- 
tal values of T c . It has therefore become widely, but not 
universally, accepted that the intramolecular H g phonons 
drive the superconductivity. Other phonon modes, such 
as librations, Cgo-Cgo and alkali-Cgo vibrations are be- 
lieved to play a small role for the superconductivity. It 
has, however, also been argued that an electronic mech- 
anism may drive the superconductivity (Chakravarty et 
al, 1991; Baskaran and Tosatti, 1991; Friedberg et al, 
1992). 

The attractive interaction between two electrons in- 
duced by the electron-phonon interaction is small (<~ 
1/10 eV) compared with the Coulomb repulsion between 
two electrons on the same Cgo molecule (U ~ 1 — 1^ 
eV). Nevertheless, for conventional superconductors it is 
argued that the effects of the Coulomb interaction are 
drastically reduced by retardation effects, due to the very 
different energy scales for phonons and electrons. The 
dimensionless Coulomb pseudopotential /z* is therefore 
believed to be drastically renormalized for such systems. 
For doped Cqo solids it has been controversial whether 
the intramolecular (Ej) or intermolecular (W) (Ander- 
son, 1991) energy scale is relevant for the retardation 
effects. In the latter case one may expect the retarda- 
tion effects to be small, since u p h <~ W. It has there- 



fore been asserted that the phonons alone cannot explain 
the superconductivity in A 3 Cgo (Anderson, 1991). Later 
work has provided support for the intermolecular energy 
scale being the relevant one, but nevertheless found that 
the electron-phonon mechanism may be sufficient if the 
Coulomb interaction is screened as efficiently as predicted 
by RPA (Gunnarsson and Zwicknagl, 1992). 

In the A3C60 (A=K, Rb) compounds many energy 
scales are similar co p h ~ Ejt ^ w p i ~ W < U, where Ejt 
is the Jahn- Teller energy of C^~ 3 , which may be of the or- 
der of several tenths of an eV. Furthermore, certain cou- 
pling constants are intermediate (A ~ 1), (g/u p i) 2 ~ 1). 
This is a situation where we may expect particularly in- 
teresting physics, but also a situation which is partic- 
ularly difficult to treat theoretically. W < U suggests 
very strong correlation effects (Lof et al., 1992), which 
have, however, not been considered in most treatment of 
fullerides. The similar magnitude of u p h and W raises 
doubt about the strength of the retardation effects, as 
discussed above, and it implies that Migdal's theorem is 
not valid. Corrections to Migdal's theorem have been 
considered in a few cases (see, e.g., Pietronero al., 1995, 
Takada 1993). Jahn- Teller- like effects, beyond what is in- 
cluded in the Migdal theory, may also play a role (Aucr- 
bach et al, 1994; Manani et al, 1994). The fact that 
Ejt may not be much smaller than W could also be im- 
portant for the insulating character of A4C60 (Fabrizio 
and Tosatti, 1996). In view of this, it is not surpris- 
ing that many of the normal state properties are unusual 
(see, e.g., Sawatzky (1995)) and several interesting, un- 
resolved questions remain. 

There have been a number of reviews of the fullcrcnes. 
Several articles appeared in J. Chcm. Phys. Solid Vol. 
53, No. 12 (Fischer and Cox, 1992), in Vol. 30 No. 8 of 
Carbon (1992), in Vol. 48 of Solid State Physics (1994), 
e.g., by Pickett (1994) and by Licber and Zhang (1994), 
and in a book edited by Billups and Ciufolini (1993). 
Early reviews were given by Weltner et al. (1989), 
by Kroto et al (1991) by Wilson et al (1992) and 
by Loktev (1992). A number of conference proceedings 
have appeared (Hammond and Kuck, 1992; Kumar et al, 
1993; Fischer et al, 1993; Prassides et al, 1994, Kuz- 
many et al, 1993; 1994a; 1995a). The field has further 
been reviewed by Holczcr and Whetten (1992) and by 
Ramirez (1994) from an experimental and by Gclfand 
(1994) from a theoretical point of view. A general re- 
view with the emphasis on the superconductivity was 
given by Hebard (1992). Specific topics have recently 
been review by Golden et al. (1995) (high energy spec- 
troscopies), Kuzmany et al. (1995b) (infrared spectro- 
scopies), Pintschovius (1996) (neutron scattering) Kuz- 
many et al (1994b) (Raman scattering) and Buntar and 
Weber (1996) (magnetic properties). The theory of su- 
perconductivity in general has been presented in several 
text-books, e.g., Schrieffer (1964), Parks (1969) and Tin- 
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kham (1975), and in reviews, e.g., Allen and Mitrovic 
(1982) and Carbotte (1990), and it will not be addressed 
here. Popular descriptions of the fullerenes have been 
given by Baggot (1994) and Aldersey- Williams (1995). 

TABLE I. Structure, valence and band involved for certain 
classes of fullerene superconductors. (AE = Alkaline Earth, 
RE = Rare Earth) 



Class 



Structure 



Valence 



A3C60 

Na 2 A 2 ,C 6 o 

As-^Ba^Ceo 

(AE^CeO 

(RE) x C6o 



Fm3m 

Pa3 

Fm3m 

Im3, sc, orthorhomb 
orthorhomb 



n = 3 (ti u ) 

n < 3 (t lu ) 

n > 3 (ti u ) 

n > 6 (tig) 



II. EXPERIMENTAL RESULTS 



A. Transition temperature 



The superconducting fullerides may be divided in five 
classes; 1) A 3 Cgo (A = K, Rb, Cs or some combination of 
these elements, 2) Na 2 A a; C 6 o (x < 1) ( A = K, Rb, Cs or 
a combination of these, 3) Aa-^Ba^Cgo (A=K, Rb, Cs), 
4) Ca^Cgo (x ~ 5) (Kortan et al, 1992a) and BagCeo 
(Kortan et al, 1992b) or Ba4Ceo (Baenitz et al, 1995) 
and SreCeo (Kortan et al, 1994) and 5) Yb2.7sC6o and 
Sm^Ceo- A summary is given in Table I. 

These classes differ in particular in their structure and 
the C 60 valence. In the first class and in the second class 
with x = 1 three electrons are doped into the t\ u band, 
while in the second class with x < 1 the doping is less 
than three electrons. In the third class the filling of the 
t\ u band is larger than three and in the fourth class the 
ti u band is full and the ti g band is populated (Chen et 
al, 1992; Wertheim et al, 1992; Erwin and Pederson, 
1993; Saito and Oshiyama, 1993; Knupfer et al, 1994). 
The tig states furthermore strongly hybridize with alka- 
line earth states (Chen et al, 1992; Erwin and Pederson, 
1993; Saito and Oshiyama, 1993), in contrast to the very 
weak hybridization of the ti u states with the alkali states 
in A 3 C 60 . 

In the first three classes, the Ceo molecules are located 
on a fee lattice. The first (Stephens et al, 1991) and third 
(Yildirim, 1996b) class have the space group Fm3m, and 
the molecules more or less randomly take one of two pref- 
erential orientations. In the second class the space group 
is Pa3 (Kniaz et al, 1993; Prassides et al., 1994), pos- 
sibly with some distortions for Na2KC6o and Na2RbCgo 
(Maniwa et al, 1995). For these systems there are four 
Ceo molecules per unit cell, all having different orien- 
tations. In the fourth group CasCeo has a simple cu- 
bic structure with the molecules at the fee lattice posi- 
tions (Kortan et al., 1992a), Ba4Ceo has an orthorhombic 



structure (Baenitz et al., 1995), BagCgo (Kortan et al, 
1992b) and Sr 6 Cgo (Kortan et al., 1994) have a bee struc- 
ture (Im3). 

The T c as a function of the lattice parameter a has 
been determined by a large number of groups for both 
the first (Tanigaki et al., 1991; Sparn et al., 1991; Flem- 
ing et al, 1991; Zhou et al, 1992) and second (Tanigaki 
et al, 1992; Yildirim et al., 1994; Mizuki et al, 1994, 
Yildirim et al., 1995) (for x = 1) class. In Fig. 3 these 
results are summarized. The figure illustrates that in 
the first class, T c varies smoothly with a, relatively in- 
dependently of the alkali atoms in the compound. Thus 
if pressure is applied to Rb 3 C6o to give the same a as 
for K3C60 at normal pressure, the corresponding T c 's are 
approximately equal (Zhou, 1992). It is therefore the 
lattice parameter, not the alkali atom, which determines 
T c . This seems to exclude the alkali vibrations as the 
main driving force for the superconductivity, since the 
vibration frequency and thereby T c would then be ex- 
pected to depend on the mass of the alkali atoms. Re- 
cent work has, however, questioned these results, finding 
a different pressure- volume relation for Rb 3 C6o than the 
earlier work and thereby a different T c — a dependence 
(Diederichs et al, 1996). A certain deviation between 
Rb2CsC6o under pressure and Rb 3 C6o was also observed 
by Schirbcr et al. (1993). This would imply that chang- 
ing the lattice parameter by physical or chemical pressure 
does not lead to identical results, although the deviations 
are not dramatic. A smooth variation with lattice param- 
eter is observed for the second group (x = 1) (see, Fig. 
3). It is, however, interesting that the slopes of the two 
classes are very different. 

(NH 3 ) 4 Na2CsC6o probably also belongs to the first 
class (with the likely structure Fm3 or Fm3m) and its 
T c — 29.6 K is essentially what one would expect from 
the large lattice parameter a = 14.473 A (Zhou et al, 
1993). Cs 3 C6o has a different (noncubic) structure and 
a different pressure dependence than the other A 3 C@o 
compounds and it should not be considered as a meme- 
ber of the first class. The maximum T c found for Cs 3 Ceo 
is 40 K (Palstra et al., 1995). Another interesting sys- 
tem is NH 3 K 3 C6o which like K 3 C6o has three electrons 
in the t iu orbital but an orthorhombic distortion of the 
fee lattice (Rosseinsky et al., 1993; Zhou et al., 1995). At 
ambient pressure this system shows no superconductivity 
above 4 K and it is possibly an insulator. Under pres- 
sure, however, superconductivity occurs with the maxi- 
mum observed T c being 28 K. It had been suggested that 
only cubic fullerides can be superconductors, but Cs 3 C6o 
and NH 3 K 3 C6o under pressure now provide examples of 
noncubic superconducting fullerides. These compounds 
nevertheless raise interesting questions about the role of 
the deviation from cubic symmetry, since the application 
of pressure reduces this symmetry lowering and increases 
T c (Palstra et al, 1995; Zhou et al, 1995). Alternatively, 
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the closeness to a Mott transition may be important for 
the pressure dependence (Palstra et al, 1995; Zhou et 
al, 1995). 
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FIG. 3. T c as a function of the lattice parameter a for the 
Fm3m (triangles, circles and open squares) and Pa3 (filled 
squares) superconductors. For K3C60 and RD3C60 the lat- 
tice parameter was varied by changing the pressure, while for 
M^M^Ceo (with M, M'=K, Rb, Cs) the lattice parameter 
was varied by changing the composition. (After Yildirim et 
al. 1994). 



stoichiometric A3C60 is an insulator (Lof et al, 1992). 
In such a case one might expect T c to increase when the 
occupancy of the t\ u orbital moves away from integer val- 
ues, in analogy with the high T c superconductors, while 
Fig. 4 shows the opposite trend. 
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FIG. 4. T c as a function of C60 valence n for the Na2Cs^C60 
(x < 1, n < 3) and Ks-xBa^Ceo (n > 3) compounds. T c is 
scaled by the T c (T™ ax ) of the end members Na2CsC6o and 
K3C60- The heavy line is a guide for the eye. Observe the 
different lattice structures. (After Yildirim et al., 1996a). 



Yildirim et al. (1996a) have studied Na 2 Cs :E C6o for 
x = 0.25, 0.50 and 0.75. The values of T c found for these 
compounds are shown in Fig. 4 as a function of the ex- 
pected valence n = 2+x. It is found that T c drops quickly 
as the valence is reduced. For n = 2.5 no superconduc- 
tivity was observed above T = 0.5 K, while the systems 
are believed to be metallic down to n=2.25 (Yildirim et 
al, 1996a). 

Yildirim et al. (1995; 1996ab) studied Rba^^Ba^Ceo 
for x=0.25, 1 and 2. Based on Raman data and theoret- 
ical considerations, they concluded that the Ceo valence 
is n = 3 + x and that Ba is essentially completely ionized. 
The results for T c are shown in Fig. 4. T c drops rapidly 
with x and for x = 1 (n = 4) and x = 2 (n = 5) no 
superconductivity was observed above T = 0.5 K. These 
compounds show a very small Pauli susceptibility, and it 
was tentatively concluded that they are weakly metallic 
(Yildirim et al, 1996b). Fig. 4 suggests a rapid drop of 
T c when the Ceo valence deviates from 3. The reduction 
for n > 3, but not for n < 3 (see Sec. VIIA.2), may be 
explainable by the reduction of the density of states with 
increasing filling, taking the merohedral disorder into ac- 
count (Yildirim et al, 1995). The reduction of T c as n is 
reduced below 3 seems not be understood. Fig. 4 further- 
more appears to be inconsistent with the conclusion that 



In the fourth group the measured transition temper- 
atures are 8.4 K (Ca.,,C 6 o), 4 K (Sr^C^) and 7 K 
(Ba^Cgo). 

Recently, superconductivity has been found in 
RE-rCeo, where RE is a rare earth atom. Ozdas et al. 
(1995) discovered that Yb2.75C6o is superconducting with 
T c = 6 K. This system has a slight orthorhombic dis- 
tortion of the Pa3 structure and the structure is Pcab. 
Based on EXAFS measurements it was concluded that 
Yb is divalent, which would imply the Cqo valence 5.5 
(Ozdas et al, 1995). Chen and Roth (1995) found super- 
conductivity in Sm^Cgo {x < 3) with T c = 8 K. 

It is also interesting to observe that Li 2 AC6o is not su- 
perconducting (Hirosawa et al, 1994). This system has 
a different structure than Na2^4C6o and AaCeo- It has 
the space group Fm3m, but with an orientational state 
of Ceo which is best modeled as a quasi-spherical unit 
(Hirosawa et al, 1994). 



B. Superconductivity gap 

The superconductivity gap A is of great interest, since 
a value of the reduced gap 2 A /T c which is substantially 
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larger than the BCS value (3.53) indicates that strong- 
coupling effects are important. For intermediate cou- 
pling, it is found that (Mitrovic et al, 1984) (see, also 
Carbotte (1990) for a discussion) 



^=3.53[l + 12.5(^) 2 ln(^)], 



(2) 



where u>i n is the logarithmic average phonon frequency. 
This formula says that the deviation from the BCS re- 
sult is small if T c /w/„ << 1. An early measurement 
of A was performed using point contact tunneling in 
STM, giving 2A/T C = 5.3 for Rb 3 C 60 (Zhang et al., 
1991b). Later STM measurements on single crystals have 
given reduced gaps 2A/T C = 5.4 (Jess et al, 1994) and 
2A/T C = 2.0 - 4.0 (Jess et al, 1996). The variations 
in 2A/T C were ascribed to variations in the local stoi- 
chiomctry at the tip position. From measurements using 
NMR the values 3.0 for K 3 C 60 and 4.1 for Rb 3 C 60 (Ty- 
cko et al., 1992), 4.3 for K 3 C 6 o (Sasaki et al., 1994) and 
3.4 ± 0.2 K 3 C 60 (Auban-Senzier et al, 1993) were de- 
duced. Recent NMR measurements were found to be in 
good agreement with the BCS value for the gap (Stenger 
et al, 1995). Using muon spin relaxation, Kicfi et al. 
(1993) obtained the value 3.6 for Rb 3 C 6 o- Recent optical 
measurements have given the values 3.44 for K 3 C 6 o and 
3.45 for Rb 3 C 60 (Degiorgi et al., 1994, Degiorgi, 1995). 
Finally, the reduced gap was measured in photoemission, 
giving the value 4.1 for Rb 3 C 6 o (Gu et al., 1994). Wc 
note that the data obtained from different experiments 
show a substantial variation. Most recent results, how- 
ever, tend to scatter around the weak-coupling limit and 
are consistent with this limit, if a reasonable uncertainty 
is assumed for these experiments. 



C. Isotope effect 

The isotope effect may provide interesting informa- 
tion about the mechanism for superconductivity. In the 
BCS theory for a system with only one type of ions 
with the mass M, the transition temperature behaves as 
T c <~ M~ a , where a — 0.5. For the A 3 C6o compounds, 
the Ceo phonons are expected to be the important ones, 
and one is therefore interested in the isotope effect when 
12 C is substituted by 13 C. For the case of complete (99%) 
substitution, it was found that a — 0.30±0.06 for K 3 C 60 
(Chen and Lieber, 1992) and a = 0.30 ±0.05 for Rb 3 C 60 
(Chen and Lieber, 1993). Interesting and contradictory 
results have been obtained for incomplete substitution. 
Chen and Lieber (1993) found that for a single peak mass 
distribution of Cqo molecules (Rb 3 ( 13 Co.55 12 Co. 45)60), 
the value of a was consistent with the result for complete 
substitution. On the other hand, they found that a two 
peak mass distribution (Rb 3 ( 13 C6o)o.5( 12 C6o)o.5) gave a 
much larger value of a (~ 0.8). A value of a which 



is larger than 0.5 is unusual, but it was also observed 
by three other groups. Ebbesen et al (1992a) obtained 
a = 1.4±0.5 for Rb 3 C 60 (33% substitution), Zakhidov et 
al (1992) obtained a = (1.2 - 1.43) ± 0.2 for K 3 C 60 and 
a = (2 - 2.25) ± 0.25 for Rb 3 C 60 (60% substitution) and 
Auban-Senzier et al (1993) obtained a — 1.45 ± 0.3 for 
82 % substitution. However, these results were observed 
for distributions of Cgo masses with essentially only one 
peak, while the a > 0.5 was obtained by Chen and Lieber 
(1993) for a two-peak distribution. Finally Ramirez et al 
(1992a) obtained a = 0.37 ± 0.05 for mass distribution 
of 15 % 12 C 60 and the rest ( 13 Co. 9 12 Co.i) 6 o- For the par- 
tially substituted compounds, the experimental situation 
thus appears to be unclear. The results are summarized 
in Fig. 5. The theoretical considerations have so far, 
however, been made for completely substituted systems, 
where a = 0.3 should be the best available experimental 
result. The theoretical treatment of partially substituted 
compounds remains a challenging problem. 
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FIG. 5. The isotope shift AT C of T c for K 3 Ceo (squares) and 
RbsCso (diamonds) as a function of the fraction p of the 
substitution (p = 1 completely 13 C). The dashed line shows 
the BCS result a — 0.5 and the full line a — 0.3. The results 
are obtained from (4) Ebbesen et al. (1992), (5) Ramirez 
ct al. (1992a), (6) Chen and Liber (1992), (7) Zakhidov ct 
al. (1992), (8) Bethune et al. (1992), (9) Chen and Lieber 
(1993), (10) Auban-Senzier et al. (1993). (After Deaven and 
Rokhsar, 1993). 



The isotope effect with respect to the alkali dopants 
have also been studied. It was found that within the ex- 
perimental uncertainty there is no isotope effect (Ebbesen 
et al 1992b; Burk et al. 1994ab). This is consistent with 
the experimental result that applying pressure to Rb 3 Ceo 
so that the lattice parameter is reduced to that of K 3 Cgo 
leads to about the same T c for both systems (Fleming 
et al. 1991). This strongly suggests that the alkali ions 
have a very weak influence on T c except for the indirect 
influence via the lattice parameter. 
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D. Other properties 

A3C60 is a type II superconductor. The lower criti- 
cal field H c i has been measured to be 130 Oe (Holczer 
et al, 1991b) for K 3 C 60 and 120 Oe (Sparn et al, 1992) 
and 60±5 Oe (Ramirez et al, 1994) for Rb 3 C 6 o. To es- 
timate the penetration depth A(0), one can now use the 
relationship with H c i A(0) (Tinkham, 1975) 



H, 



<Po 1 , A(0) . 
47rA(0)2 m ^GL' 



(3) 



where (f>o — hc/e and £gl is the Ginzburg-Landau co- 
herence length. Putting ^gl — 30 A, which is reason- 
able in view of the results below, one obtains the val- 
ues A(0) = 2500 A for K 3 C 60 and 2600 and 3700 A for 
Rb 3 C6o- The penetration depth has also been estimated 
from /xSR with the results 4800 A for K 3 C 60 and 4200 A 
for Rb 3 C6o (Ucmura et at, 1991; 1993), and from the op- 
tical conductivity with the result 8000 A for both K 3 Ceo 
and Rb 3 C 60 (Degiorgi et al, 1992) and from NMR data 
giving 6000 A for K 3 C 60 and 4600 A for Rb 3 C 60 (Tycko et 
al, 1992). Uemura et al. (1991) also deduced the tem- 
perature dependence of penetration depth A(T). They 
could fit these results with the formula [1 — (T/T c ) a ], 
with a — 3.2. From this they concluded that K 3 C6o 
probably is an isotropic s-wave superconductor, as ex- 
pected. Based on values for A(0) and £, Uemura et al. 
(1991) also estimated an "effective" Fermi temperature 
T^ff . Combining their small value (470 K) for Tp^ with 
T c in a T c vs. Tp* plot, they found that the A 3 C6o be- 
longs to a group of "exotic" superconductors, including 
the High T c cuprates and the organic BEDT systems. 

The upper critical field H C 2(0) is of interest, since it al- 
lows an estimate of the coherence length via (Tinkham, 
1975) 



H c2 (0) 



2tt£_ G l 



(4) 



Several groups have measured H c2 = dH c ijdT and 
then obtained the value of H C 2(0) from the Wcrthcimer- 
Helfand-Hohenberg (Wertheimer et al., 1966) theory. 
Based on this, £gl was deduced to be 26 A for K 3 C6o 
(Holczer et al., 1991b) and 20 A for Rb 3 Ceo (Sparn et al, 
1992) from dc-magnctization. The high-field susceptibil- 
ity gave the values 29 A to 33 A (Boebinger et al, 1992) 
and 35 A (Holczer and Whetten, 1992) for K 3 C 60 and 
30 A (Holczer and Whetten, 1992) for Rb 3 C 60 . From 
the ratio of £gl/A it is clear that A 3 C^ is a type II 
superconductor. The magnetic properties of the Ceo su- 
perconductors have been reviewed by Buntar and Weber 
(1996). 

The jump in the specific heat at T c was estimated by 
Ramirez et al. (1992) to be AC = 68 ± 13 mJ/mole-K 2 
for K 3 C6o and by Burkhart and Meingast (1996) to be 



64±14 and 75±14 mJ/molc-K 2 for K 3 C 60 and Rb 3 C 60 , 
respectively. The Hebel-Slichter peak was not seen in 
early NMR measurements (Tycko et al., (1992), but has 
been seen in /iSR (Kiefl et al., (1993)) and more recent 
NMR (Stenger et al., 1995) measurements. The vertex 
glass and vertex fluid phase diagram of A 3 C6o has been 
studied by Lin et al. (1994). 

In Table II we have summarized some superconducting 
parameters. The table has mainly been extracted from 
Ramirez (1994), but it has been partly updated in the 
view of newer experiments. 

TABLE II. Superconducting parameters for K3C60 and 
Rb 3 C 6 o- (Mainly after Ramirez (1994)). 



Property 



K 3 C e 



Rb 3 C 6 



A {nSR) 
A (optical) 
A (NMR) 

£(> (Xac) 

AC/T C 

2A/k b T c (NMR) 
2A/k b T c (optical) 
2A/k b T c (/iSR) 



19.5 K 

4800 ± 200 A 
8000 ± 500 A 
6000 A 
29-35 A 



29.5 K 
4200 A 
8000 ± 500 A 
4600 A 
30 A 



64 ± 14 mJ/mole-K 2 75 ± 14 mJ/mole-K 2 
3.0 4.1 
3.44 3.45 

3.6 ±0.3 



2A/k b T c (tunneling) 5.3 ± 0.2 
2A/k b T c (tunneling) 



5.2 ±0.3 
5.4,2.0-4.0 



III. PHONONS AND ELECTRON-PHONON COUPLING 



The phonons and the electron-phonon interaction are 
very important for the properties of alkali-fullcrides. The 
electron-phonon interaction is usually assumed to cause 
the superconductivity. It should, however, also be impor- 
tant for the transport properties in general, and may also 
play a substantial role for other electronic properties. 

The phonons of alkali-fullerides can be divided in sub- 
groups, which reflect the molecular nature of the solid. 
The highest modes are due to intramolecular vibrations, 
with energies in the range 273-1575 cm -1 (0.034-0.195 
eV) for an undoped Cqq solid (Bcthune et al., 1991). 
The intermolecular phonons occur at substantially lower 
energies in the range up to about 140 cm -1 (17 meV) 
(Prassides et al., 1991; Pintschovius, 1996). These modes 
contain both Ceo-Ceo and A-C50 (A=K, Rb) vibrations. 
Finally, there are librations, which occur at an energy 
of about 4-5 meV (Christides et al., 1992; Schober et 
al., 1994). These various phonons are schematically il- 
lustrated in Fig. 6. It is believed that the intramolecular 
modes play the main role for the superconductivity. We 
therefore focus on these modes, but also the other modes 
are discussed. 
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FREQUENCY (cm \> 

FIG. 6. Schematic representation of various phonons in 
A3C60 compounds. The figure shows (a) librations, (b) optic 
C60-C60 phonons, (c) A-C60 phonons and (d)-(e) intramolec- 
ular H 9 modes. The figure indicates the radial and tangen- 
tial character of the low-lying and and high-lying H 9 modes, 
respectively. (After Hebard, 1992). 



A. Intramolecular modes 

For the intramolecular phonons, the effects of the in- 
teractions between the Cgo molecules are small. These 
phonons can therefore be approximately classified by the 
icosahedral point group. In alkali-fullerides, the t\ u band 
is (partly) populated. The coupling of the phonons to 
this band is therefore of particular interest. From sym- 
metry arguments, it follows that only the two phonons 
with A g and the eight phonons with H g symmetry can 
couple (Varma et al, 1991; Lannoo et al, 1991). Wc 
therefore focus our interest on these modes. 



1. Calculations of phonon frequencies 

There have been a number of calculations of the 
phonon frequencies. In Table III we show the results 
of calculations using the semi-empirical MNDO method 
(Varma et al., 1991), an empirical force-constant model 
(Jishi et al., 1992) and the ab initio density functional 
method (Hohenberg and Kohn, 1964; Kohn and Sham, 
1965) in the local density approximation (LDA) but us- 
ing different numerical methods ( Antropov et al, 1993b; 
Quong et al., 1993; Faulhaber et al, 1993). These calcu- 
lations were all for a free Cgo molecule, except the one by 
Antropov et al., (1993b), which was for K3C60 in the so- 
called uni-directional structure, where all Ceo molecules 
have the same orientation (where a cubic axis and a two- 
fold molecular axis coincide). These calculations can be 
considered as representative for the large number of cal- 
culations that have been performed (see also, e.g., Negri 
et al, 1988; Adams et al, 1991; Kohanoff et al, 1992; 
Jones et al, 1992; Onida and Benedek, 1992; Bohnen et 



al, 1995). These results are compared with experimental 
results obtained from Raman scattering (Bcthunc et al, 
1991). We can see that the best empirical and ab initio 
calculations have errors of the order one per cent. The 
small errors in, for instance, the calculations of Quong et 
al, (1993) suggests that the LDA errors for the frequen- 
cies are small. 

Bohnen et al. (1995) have calculated the changes in the 
phonon frequencies when going from solid Ceo to K3C60, 
and found a softening of the modes, in particular the 
high-frequency ones. 



TABLE III. Phonon frequencies (in cm -1 ) for the H g and A s 
modes. Different theoretical calculations are compared with 
experimental results for undoped solid Ceo- 



Mode 


Varma a 


Jishi b 


Energy 
Antropov Quong d 


Faulhaber Exp. f 


H 9 (8) 


1721 


1575 


1462 


1576 


1567 


1575 


H B (7) 


1596 


1401 


1387 


1443 


1425 


1428 


H 9 (6) 


1406 


1217 


1290 


1244 


1200 


1250 


H 9 (5) 


1260 


1102 


1091 


1093 


1067 


1099 


H 9 (4) 


924 


788 


785 


767 


750 


774 


H 9 (3) 


711 


708 


753 


727 


640 


710 


H 9 (2) 


447 


439 


454 


439 


421 


437 


H 9 (l) 


263 


269 


281 


258 


249 


273 


A 9 (2) 




1468 


1463 


1499 


1493 


1470 


A 9 (l) 




492 


458 


478 


459 


496 



a Varma et al, 1991 
b Jishi et al, 1992 
Antropov et al, 1993 
d Quong et al, 1993 
°Faulhaber et al, 1993 
f Bethune et al, 1991 

The four highest modes are primarily made up of vi- 
brations tangcntially to the Cqo surface, while the three 
lowest modes have mainly radial character (Schluter et 
al, 1992ab; Antropov et al, 1993b), as indicated in Fig. 
6. The reason is that the molecule is much stiffer against 
tangential than radial distortions, since the former re- 
quire larger bond length changes. The phonon frequen- 
cies are high, extending up to about 0.2 eV. This is due 
to both the light mass of carbon and to the stiffness of 
the C-C bond. 



2. Calculations of electron-phonon coupling 

The electron-phonon coupling plays a crucial role in 
the theory of superconductivity. We introduce the di- 
mensionless coupling to the phonon mode v (Rainer, 
1986) 

2 1 

^ = ^J2 — E i^ k ,-(k +q) H 2 

q n,m,k 
X<5(£„k)<5(£m(k+q) - £nk - Uvq), (5) 
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where iV(O) is the electron density of states per spin 
at the Fermi energy and ov q is the energy of the vth 
phonon with the wave vector q. The energy of the nth 
electronic band with the wave vector k is e„k- Finally, 
5nk,m(k+q) ( u ) is the matrix element between the two elec- 
tronic state nk and m(k + q) of the potential created 
when a phonon v<\ is excited. Thus it follows 

5nk,m(k+q)(^) = (nk| &V vq \m(\t + q)> — ;=== , (6) 

where AV„k is the induced potential per unit displace- 
ment along the phonon normal coordinate, y / 2Mw 1/q is 
the phonon amplitude and M is the reduced mass. 

The calculation of A requires a double sum over a Bril- 
louin zone of a quantity (g) which is very time consum- 
ing to calculate. In general, an accurate calculation of 
A is therefore very complicated. For the intramolecu- 
lar modes in the fullerides, important simplifications fol- 
low from the large difference between the intramolecular 
(Ej) and intermolecular (W) energy scales (Lannoo et 
al, 1991). One can then assume that the phonon eigen- 
values and eigenvectors are independent of q. If in ad- 
dition we assume that W << Ej in the calculation of 
5nk.m(k+q) i v ) , but keep the full structure of |nk >, the 
complicated Brillouin zone sums can be performed an- 
alytically, and the problem is essentially reduced to a 
molecular problem (Lannoo et al., 1991). The accuracy 
of these assumptions has been tested for a simple tight- 
binding model, where it was possible to perform the Bril- 
louin zone summations (Antropov et al., 1993b). It was 
found that the simplified calculation agreed fairly well 
with the full calculation, justifying the assumptions of 
Lannoo et al. (1991). 

Thus it is sufficient to calculate the shift Ae ua of the 
t\ u level a for a free Cgo molecule per unit displacement 
of the vth phonon coordinate. One then finds that 

va v 

This gives a molecular specific quantity which must then 
be multiplied by the density of states iV(0), which is de- 
termined by the intermolecular hopping. The molecular 
calculation therefore gives a value for X l/ /N(Q), and it is 
appropriate to compare this quantity from different cal- 
culations. We observe that the quantity X v /N(0) grows 
with the intramolecular energy scale, while iV(0) is in- 
versely proportional to the intermolecular hopping. The 
large difference between these two energy scales there- 
fore favors a large value of A (Schluter et al, 1992a). It 
is also clear that as the lattice parameter increases, the 
overlap between the Ceo molecules is decreased, reducing 
the band width. This leads to a larger density of states 
-/V(O) and a larger A, explaining the increase in T c with 
lattice parameter (Rosseinsky et al., 1991; Flemming et 



al, 1991; Tanigaki et al., 1993; Ramirez, 1994) (see, also 
Fig. 3). 

We notice that the results are quite sensitive to the 
phonon eigenvectors, and small changes in the eigenvec- 
tors can lead to large shifts of the coupling strength be- 
tween different phonons. One can show, however, that 
^^(Ae^) 2 is independent of the phonon eigenvectors 
(Antropov et al., 1993b). A shift of coupling strength 
between phonons with similar energies therefore has a 
small effect on the total value of A, while a shift be- 
tween phonons of very different energies can dramatically 
change the total A due to the energy denominator in Eq. 
(7). For instance, in the calculation of Antropov et al. 
(1993b), mixing 5 % of the weight (VoM) of the sev- 
enth eigenvector into the eighth and vice versa changes 
the couplings 0.020 and 0.022 in Table IV to 0.010 and 

0. 030, i.e. a small change in the total coupling. The same 
mixing between the second and the eighth eigenvectors 
changes the couplings 0.006 and 0.022 to 0.033 and 0.019, 

1. e. a large change of the total coupling. 

TABLE IV. The partial electron-phonon coupling constants 
A„/V(0) (in eV) according to different theoretical calcula- 
tions. The energies (in cm -1 ) of the modes for the undoped 
system are shown. 



Mode Energy A„/iV(0) 

Varma a Schluter b Antropov c Faulhaber d 



H 9 (8) 


1575 


.011 


.009 


.022 


.009 


H 9 (7) 


1428 


.034 


.013 


.020 


.015 


H 9 (6) 


1250 


.000 


.003 


.008 


.002 


H 9 (5) 


1099 


.006 


.001 


.003 


.002 


H 9 (4) 


774 


.000 


.007 


.003 


.010 


H 9 (3) 


710 


.001 


.004 


.003 


.001 


H 9 (2) 


437 


.001 


.007 


.006 


.010 


H 9 (l) 


273 


.003 


.008 


.003 


.001 


E h 9 




.056 


.052 


.068 


.049 



a Varma et al, 1991 
b Schluter et al, 1992b 
c Antropov et al, 1993b 
d Faulhaber et al., 1993 



In Table IV we show some typical results of calcula- 
tions of X/N(0). The calculation of Varma et al. (1991) 
was based on the semi-empirical MNDO method. The 
calculation of Schluter et al. (1992b) used an empirical 
bond charge model for the phonons and a LDA calcu- 
lation for the electronic part. The work of Antropov et 
al. (1993b) and Faulhaber et al. (1993) used ab ini- 
tio LDA calculations both to generate the phonons and 
to calculate their interactions with the electrons. There 
are several other calculations. Mazin et al. (1992) per- 
formed an ab initio LDA calculation in the ASA-LMTO 
formalism and used an empirical phonon model to obtain 
A/V(0) = 0.037 cV. Coulon et al. (1992) performed an 
ab initio LDA calculation using a pseudopotential Gaus- 
sian method and found X/N(0) > 0.04 eV. Asai and 
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Kawaguchi (1992) obtained X/N(0) = 0.070 eV. The cou- 
pling to the A g (2) mode for a free molecule was estimated 
by Stollhoff (1991) to be X Ag{2) /N(0) = 0.0025 eV and by 
Scherrer and Stollhoff (1993) to be X Ag{2 )/N(0) = 0.01 
eV. Relatively small values for the coupling to the A s (2) 
mode were also deduced by Schliiter et al. (1992a) 
(^A g (2)/N(0) = 0.005 eV), Antropov et al. (1993b) 
{X Ag{2 )/N(0) = 0.009 eV) and Pickett et al. (1994) 
{X Ag{2 )/N(0) = 0.002 eV). These calculations do not in- 
clude the screening effects discussed below. The clcctron- 
phonon interaction has also been calculated by Jishi and 
Dresselhaus (1992) and by Chen and Goddard (1993). 

Table IV illustrates that there are substantial differ- 
ences between the different calculations in terms of the 
distribution of coupling strength between the different 
modes. Thus Varma et al. (1991) found that almost 
all the coupling strength is to the two highest modes. 
Antropov et al. (1993b) found more coupling strength 
to the lower phonons and this was still more true for the 
calculations of Schlutcr et al. (1992b) and Faulhaber et 
al. (1993). Nevertheless, all four calculations find the 
strongest coupling to one of the uppermost phonons. The 
large variation in the distribution of coupling between the 
different phonons in Table IV may to a substantial ex- 
tent be due to the sensitivity to the phonon eigenvectors 
discussed below Eq. (7). These variations suggest that 
there are substantial uncertainties in the calculations. 

In estimates of the total A, the contributions from the 
A g modes are usually not included. The reason is that 
the A g contribution is expected to be screened out in 
A 3 C 60 (A=K, Rb) (Schluter et al, 1992b). This is illus- 
trated by Raman scattering results for A3C60 (A=K, Rb) 
(Duclos et al, 1991, Pichlcr et al, 1992). As discussed 
in Sec. III. 3, the additional broadening of the A g modes 
in A 3 C 6 o compared with Ceo should be a measure of the 
electron-phonon coupling. Experiments, however, show 
little or no extra broadening of the A ff (2) in the metallic 
phase. For instance, Pichler et al. (1992) estimated the 
additional broadening to be only about 3 cm -1 , implying 
Ayig(2)/A(0) ~ 0.0002, i.e., an extremely small coupling. 

The corresponding additional broadening for many of 
the Hg modes is, on the other hand, large, indicating a 
much stronger coupling for these modes. The symmet- 
ric Ag mode leads to the same displacement of all the 
t\ u levels on a given molecule, and this displacement can 
largely be compensated by a transfer of charge to this 
molecule. This screening is much less efficient for the 
Hg modes, since these do not shift the center of gravity 
of the ti u levels. Simple estimates within RPA and as- 
suming static screening, suggest that the screening could 
reduce the coupling to the A g mode by one to two or- 
ders of magnitude (Antropov et al, 1993b). It would be 
interesting to investigate nonadiabatic effects, going be- 
yond static screening, since the A g (2) phonon energy is 
comparable to the energy scale of the screening. 



TABLE V. The density of states N(0) (per eV and spin) 
for K3C60 and RD3C60 according to different band structure 
calculations. Some of the results have been interpolated by 
Gelfand (1994). 



N(0) (K 3 C 6 o) N(0) (Rb 3 C 6 o)/iV(0) (K 3 C 6 o) 

6.6 a 1.26 a 
7.0 b 

8.6 C 1.18° 

8.6 d 1.21 d 

9.0 C 1.14° 

9.8 f 1.27 f 



a Erwin and Pickett (1991; 1992) 

b Antropov (1992) 

c Huang et al. (1992) 

d Satpathy et al. (1992) 

°Novikov et al. (1992) 

f Troullier and Martins (1992) 

To obtain values of A, we also need estimates of the den- 
sity of states at the Fermi energy N(Q). There have been 
several theoretical calculations of A^(0) for K 3 C6o- Full 
potential, linear combination of Gaussian orbitals gave 
iV(0) = 6.6 states/eV-spin for the socalled uni-directional 
structure, where all the C@o molecules have the same ori- 
entation and there is one Cqo molecule per unit cell 
(Erwin and Pickett, 1991; 1992). An ASA-LMTO calcu- 
lation gave N(0) = 8.6 states/eV-spin (Satpathy et al, 
1992) and a full potential LMTO gave 7.0 state/eV-spin 
(Antropov, 1992). These and some other results are col- 
lected in Table V. A calculation for the bi-directional 
structure, where the molecules have two different direc- 
tions, gave N(0) = 7.5 states/eV-spin (Antropov, 1992). 
The most relevant calculation is probably, however, for 
the orientationally disordered structure. This has only 
been performed within a tight-binding approach (Gelfand 
and Lu, 1992ab; 1993). All the structures in N(e) are 
then smeared out, and N(0) is slightly larger than 3/W, 
where W is the band width. If we assume W <~ 0.6 
eV, as many band structure calculations find (Erwin and 
Pickett, 1991; 1992; Antropov, 1992), this gives AT(0) of 
the order 5-6 states/eV-spin. 

Experimentally, the density of states has been esti- 
mated using a number of different techniques. Many es- 
timates are based on the NMR relaxation rate. It has 
often been assumed that the Fermi contact interaction 
dominates the relaxation, and large values of N(0) were 
then obtained, ft has, however, been shown that the 
relaxation is dominated by the spin-dipolar interaction, 
due to the p z character (<~ 95 %) of the states close to 
the Fermi energy (Antropov et al., 1993a). These re- 
sults were supported by Tycko et al. (1993), who found 
that the spin-dipolar coupling dominates for RbCgo and 
CsC 60 . Antropov et al. (1993) deduced that N(0) = 7.2 
and 8.1 states/eV-spin for K 3 Cgo and Rb 3 C6o, respec- 
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tively. These values are used in the following. 

N(0) has also been estimated from the specific heat 
jump AC at the superconductivity transition (Ramirez 
et al., 1992a; Burkhart and Meingast, 1996). Using the 
formula 

AC 

— ~V(0)(l + A)(1.43 + O(A 2 )), (8) 

J c 

and A = 0.068iV(0) (Antropov et al., 1993b) Burkhart 
and Meingast (1996) obtained N(0) = 6.6 and 7.4 
states/eV-spin for K 3 Ceo and Rb 3 C6o, respectively. Us- 
ing instead A = 0.147 AT(0) (Gunnarsson et al, 1995), 
they found N(0) = 5.3 and 5.9 states/eV-spin for K 3 C 6 o 
and Rb 3 C@o, respectively. 

The density of states has furthermore been estimated 
from Pauli contribution to the susceptibility obtained 
from SQUID and EPR measurements, as in Table VI. 
The SQUID results were corrected for the diamagnetic 
contribution from the cores. The result of Wong et 
al. (1992) was also corrected for the Landau orbital 
diamagnetism, while Ramirez et al. (1992b) argued 
that this contribution is negligible in the present sys- 
tems. The EPR measurements contain no diamagnetic 
contributions. The susceptibility may have an apprecia- 
ble Stoner enhancement (Ramirez et al., 1992b). Based 
on LDA calculations, this enhancement factor was esti- 
mated to be 1.34 and 1.42 for K 3 C6o and Rb 3 C6o, re- 
spectively (Antropov et al., 1993a). The accuracy of the 
LDA in this context is, however, not known. This Stoner 
enhancement has not been taken into account in Table 
VI, and the numbers in the table should be divided by 
the appropriate Stoner enhancement. The susceptibility 
has recently been calculated for a three-fold degenerate 
Hubbard model, describing the t\ u band and the on-site 
Coulomb interaction U. A quantum Monte Carlo lattice 
method was used to obtain the enhancement of the Pauli 
susceptibility due to the Coulomb interaction (Aryaseti- 
awan et al., 1996). For typical values of U ~ 1.2 — 1.3 eV, 
an enhancement of the order of a factor three was obtain. 
This would reduce the values of N(0) estimated from the 
susceptibility and roughly bring them in line with other 
estimates. 

We observe that it is hard to estimate N(0) from pho- 
toemission or yield spectroscopy (De Seta and Evange- 
listi, 1995), due to a large broadening of the band, which 
is probably caused by phonon and plasmon satellites 
(Knupfer et al., 1993) or other many-body effects. Vari- 
ous experimental estimates of N(0) have been discussed 
by Gelfand (1994). 

Table VII shows the contribution from the H 3 modes 
to A using the values of N(0) deduced from NMR, 
i.e., N(0) = 7.2 and 8.1 states/eV-spin for K 3 C@o and 
Rb 3 C6o, respectively. 



TABLE VI. The density of states JV(0) (per eV and spin) 
for K3C60 and RD3C60 as deduced from susceptibility mea- 
surements. The results have not been corrected for the Stoner 
enhancement, which would lead to reduced estimates of N(0). 

N(0) (K 3 C 6 o) N(0) (Rb 3 C 6 o)/iV(0) (KgCgo) Method 



14 a 1.36 a " " SQUID 

16 b SQUID 

ll b EPR 

15 c 1.40 d EPR 

10 d EPR 



a Ramirez et al. (1992b) 
b Wong et al. (1992) 
Tanigaki et al. (1995) 
d Wang et al. (1995) 



TABLE VII. The H 9 contribution to A according to different 
theoretical calculations, using N(0) = 7.2 and 8.1 states per 
eV and spin for K3C60 and RbsCeo, respectively. 



Mode 






A„ 






Varma a 


Schluter b 


Antropov 


Faulhaber d 


K3C60 


0.40 


0.37 


0.49 


0.35 


RbsCeo 


0.45 


0.42 


0.55 


0.40 



a Varma et al., 1991 
b Schluter et al, 1992 
c Antropov et al, 1993 
d Faulhaber et al., 1993 



It is interesting to compare the alkali-fullcridcs with 
intercalated graphite, since the latter have much lower 
values (one to two orders of magnitude) of T c than the ful- 
lerides. It has been reported that calculations for inter- 
calated graphite gives a substantially lower value (factor 
five) of A than for the fullcrides (Schluter et al., 1992a). 
This is not ascribed to a larger density of states in the 
fullcrides but to the curvature of the Ceo molecule (Ben- 
ning et al., 1991; Schluter et al., 1992ab). In intercalated 
graphite, to first order there is no coupling to displace- 
ments perpendicular to the graphite planes. In Cgo, on 
the other hand, the curvature mixes in 2s character in the 
primarily radial 2p orbital pointing out from the molecule 
and dominating the states in the t\ u band. This allows 
for a coupling of these states to radial (transverse) modes. 
Alternatively, we can think of the displacements of the 
nuclei as setting up an approximately p-like perturbing 
potential close to the nuclei. Such a potential can cou- 
ple a 2s and a 2p orbital but not two 2p orbitals on the 
same site. Pietronero and Strassler (1992) suggested as 
an alternative explanation that the violation of Migdal's 
theorem in A 3 Ceo may cause the large T c in these com- 
pounds. 
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3. Raman and neutron scattering 

Due to the electron-phonon coupling, a phonon can de- 
cay in an electron-hole pair of the appropriate energy. In 
the metallic fullerides A3C60 this decay channel is avail- 
able and contributes to the width of the phonons. This 
width is therefore a measure of the electron-phonon cou- 
pling (Allen, 1972; Allen, 1974) and for the vth phonon 
it is given by 

7 „ = 2 ^ ( ° )A -, (9) 
9v 

where 7^ is the full width at half maximum (FWHM), ui v 
is the phonon frequency, N(0) is the density of states per 
spin, A„ the electron-phonon coupling and g v is the de- 
generacy of the phonon (5 for H s ). Due to the factor cj^, 
the width of the high-lying phonons is often very large 
and those phonons are hard to observe experimentally. 

The width can be measured using Raman or neutron 
scattering. In an ordered system, Raman scattering sam- 
ples the decay of a phonon for q w and it may be mis- 
leading for the decay in the rest of the Brillouin zone. 
In A3C60 (A=K, Rb) there is, however, a strong ori- 
entational disorder (Stephens et al., 1991) and Raman 
scattering should be relevant for the determination of 
the electron-phonon coupling integrated over q (Schluter, 
1992c). The observed widths may therefore include some 
dispersion effects. These are, however, expected to be 
small for the intramolecular modes, due to the weak cou- 
pling between the molecules (Belosludov and Shpakov, 
1992). 

The widths of some of the H s phonons were measured 
very early by Mitch et al. (1992a, 1992b). They per- 
formed Raman measurements for ultrathin films of alkali- 
doped A^Ceo compounds. Due to the small thickness of 
the film, they were able to continuously dope the system 
with Na, K, Rb, Cs. For the H s (2) mode they observed 
a half-width of 14, 22, 22 and 20 cm" 1 for A=Na x=l, 
A=K £=2.1, A=Rb £=1.9 and A=Cs x—2, respectively. 
They assigned the rather ^-independent broadening in 
NaCgo to disorder, and assumed that the disorder is sim- 
ilar in the other compounds. The additional broaden- 
ing in these compounds was assigned to the electron- 
phonon coupling. Assuming that the two mechanisms 
add quadratically, they deduced the widths 17 cm^ 1 for 
K and Rb doped Cgo- To obtain a value for A2 from Eq. 
(9), one further has to assume a value for N(0). Since 
different groups have used different values for N(0), we 
have here uniformly used N(0) = 7.2 eV -1 per spin, de- 
duced from NMR for K3C60 (Antropov et al, 1993a), to 
be able to compare the different estimates. The values 
quoted below therefore differs from the values of A given 
in the original papers. In a few cases below we have also 
corrected for an incorrect prefactor in Eq. (9) which has 



appeared in the literature. In Table VIII we show differ- 
ent experimental estimates of A„/iV(0) for K 3 C6o- 

The phonon width has also been studied in neutron 
scattering for Ceo, K3C60 and RbgCgo by Prassides et al. 
(1992). By comparing the broadening for the metallic 
K3C60 and the insulating Ceo, the additional broadening 
due to the electron-phonon coupling was deduced. This 
was performed for the four lowest H s modes. The result 
for H s (2) is in good agreement with the result of Mitch 
et al. (1992b). 

Recently, the Raman spectrum was studied by Win- 
ter and Kuzmany (1996) for single crystals and at low 
temperatures. They found that the peaks for several of 
the H 9 phonons were strongly asymmetric and could be 
fitted with five components. This is shown in Fig. 7. 

TABLE VIII. The partial electron-phonon coupling con- 
stants X v /N(0) [in eV] according to different experimental 
estimates. We also show the energies (in cm -1 ) of the modes 
for the undoped system. 



Mode 


Energy 


Raman a 


A„/W(0) 
Neutron b 


Raman c 


PES d 


H 9 (8) 


1575 






.002 


.023 


H 9 (7) 


1428 






.003 


.017 


H s (6) 


1250 






.001 


.005 


H 9 (5) 


1099 






.001 


.012 


H s (4) 


774 




.005 


.002 


.018 


H 9 (3) 


710 




.001 


.002 


.013 


H 9 (2) 


437 


.022 


.023 


.014 


.040 


H 9 (l) 


273 




.014 


.034 


.019 



a Mitch et al, 1992b 
b Prassides et al, 1992 
c Winter et al, 1995 
d Gunnarsson et al, 1995 



It is interesting that the component with the largest shift 
relative to the position in the undoped crystal also has 
the largest width. This suggests that the different com- 
ponents have different electron-phonon couplings. The 
origin of such a difference in the coupling, if it exists, 
is however not clear. It is also interesting that the two 
highest H g modes seemed to have a different pattern of 
splitting. The observed line shape therefore poses a very 
interesting theoretical problem. If it is assumed that the 
formula of Allen (1972, 1974) is still applicable to each 
component by itself, the results in Table VIII are ob- 
tained. 

The measurement of the phonon broadening provides 
rather direct information about the electron-phonon cou- 
pling. We note, however, that the interpretation (Allen, 
1972) is based on Migdal's theorem. Although this 
should be an excellent approximation for most systems, 
it is questionable for the fullerides since the phonon 
frequencies and the t\ u band width are comparable 
(Pietronero and Strassler, 1992; 1994; Pietronero et al, 
1995; Grimaldi et al., 1995ab). There are also questions 
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if Jahn- Teller like effects studied for free molecules (Auer- 
bach et al., 1994; Manini et al, 1994) and going beyond 
Migdal's theorem, could survive in the solid and influence 
the results. Such higher order effects are very important 
for the photoemission spectrum of a free molecule dis- 
cussed below. In particular, the interesting line shape 
obtained by Winter and Kuzmany (1996) raises questions 
about our understanding of how to analyze the spectra. 

Reviews of Raman and neutron scattering studies of 
the fullcrides have been given by Kuzmany et al. (1994b) 
and by Pintschovius (1966), respectively. 
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FIG. 7. Raman spectrum of the H 9 (2) mode excited with a 
green laser (a) and with a red laser (b). The full lines show 
a fit to five Voigtian profiles (after Winter and Kuzmany, 
1996). 

4. Photoemission 

An alternative way of deducing the electron-phonon 
coupling is the use of photoemission data. Because of 
the relatively strong electron-phonon coupling, we ex- 
pect to see satellites due to the excitation of phonons. 
The weight of the satellites should then give information 



about the strength of the coupling. This is essentially 
the Franck-Condon effect, except that it is important to 
take the Jahn- Teller effect into account for Cgo- The 
photoemmision spectra of K3C60 and RbjCgo have been 
analyzed along these lines (Knupfer et al, 1993). Due 
to the broadening effects in a solid and due to the com- 
plications in the theoretical treatment of bands with dis- 
persion, it was not possible, however, to derive reliable, 
quantitative values for the electron-phonon coupling. 

Photoemission experiments have, however, also been 
performed for free C^~ molecules, where the theoretical 
treatment is substantially simpler (Gunnarsson et al, 
1995). In these experiments a beam of Cg ions was 
created and a photoemission experiment was performed 
using a laser light source (flu) =4.025 eV) and a time 
of flight spectrometer. To analyse the results, a model 
was used which includes the t\ u level from which the 
electron is emitted, the two A 9 and eight, five-fold de- 
generate H g modes and their coupling to the electronic 
state. For this model the ground-state can be calculated 
by numerical diagonalization to any desired accuracy. 
Within the sudden approximation (Hedin and Lundqvist, 
1969), the photoemission spectrum can furthermore eas- 
ily be calculated. A set of coupling constants are then 
assumed and the resulting spectrum is compared with ex- 
periment. The coupling parameters are then varied un- 
til good agreement with experiment is obtained, thereby 
providing an estimate of the couplings. The resulting 
spectrum is compared with experiment in Fig. 8 and the 
corresponding parameters are shown in Table VIII. 

This approach contains mainly three uncertainties. 
Firstly, it is based on the sudden approximation. This 
approximation should become very accurate as the ki- 
netic energy of the emitted electron becomes large. Due 
to the small photon energy used, to obtain a good reso- 
lution, the accuracy of the sudden approximation is un- 
clear. Spectra taken at different photon energies did not, 
however, show any systematic variation, suggesting that 
the sudden approximation may not be a serious problem 
here. We further observe that the weight of the satel- 
lites must go to zero as the photon energy is reduced 
sufficiently. If this happens in a monotonous way, the 
use of the sudden approximation would underestimate 
the couplings. As shown in Table VIII, the couplings 
derived from photoemission, however, tend to be larger 
than other estimates. The second uncertainty is more 
technical, and is related to the fact that for the resolu- 
tion available, it is not possible to distinguish between the 
coupling to the A g modes and the H s modes with similar 
energies. The couplings to the A g modes were therefore 
taken from a calculation (Antropov et al., 1993b). With 
this assumption the couplings to the H g modes can then 
be determined uniquely. An equally good fit to exper- 
iment could, however, be obtained if, for instance, the 
coupling to the A s (2), \a 9 (2)/N(0), was allowed to vary 
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between 0.00 and 0.03 eV, provided that the total cou- 
pling to the H g (7) and H s (8) modes were simultaneously 
varied between 0.07 and 0.00 eV. Thirdly, the coupling is 
determined for C^~ , while we are interested in CgQ~ in a 
solid. Although the additional charging of Cgo is not ex- 
pected to have large effects, due to the sensitivity of the 
coupling to the phonon eigenvectors, the coupling could 
still be influenced. 
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FIG. 8. The experimental (dots) and theoretical (full line) 
photoemission spectrum of Cg . The theoretical no loss 
(dashed), single loss (dotted) and double loss (dashed- 
dotted) curves are also shown. The contributions of the 
different modes to the single loss curve are given by bars 
(H 9 : open, A 9 : solid). The inset shows the experimental 
spectrum over a larger energy range (After Gunnarsson et 
al, 1995). 



5. Resistivity 



The electric resistivity can be caused by various scat- 
tering mechanisms. An important mechanism is often 
the electron-phonon interaction, where electrons are scat- 
tered under the simultaneous excitation of phonons. An- 
other mechanism is provided by the electron-electron 
scattering. The contribution from both these mecha- 
nisms goes to zero as T — > 0. Deviations from periodicity, 
e.g., the orientational disorder, on the other hand, leads 
to an essentially temperature independent resistivity. If 
the dominating contribution to the T-dcpcndcnt part of 
the resistivity comes from the electron-phonon scatter- 
ing, the resistivity provides an interesting check of the 
values of the electron-phonon coupling constants. 



O 
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CM 




FIG. 9. The experimental (dots)(Xiang et al, 1992) and the- 
oretics 1 (full curve) T-dependent resistivity p(T) normalized 
to the value at T = 260 K. 

In Fig. 9 typical results for the temperature de- 
pendence are shown. For low T, it is found that the 
temperature-dependent part of p(T) at constant pressure 
is approximately quadratic in T (Crespi, 1992). This sug- 
gests that the mechanism at these values of T is electron- 
electron scattering (Palstra et al, 1994). To obtain a rea- 
sonable magnitude, it is, however, necessary to assume 
that the matrix element for electron-electron scattering is 
of the order 1 eV, i.e., comparable to the t\ u band width. 
Such a large matrix element would suggest a rather ineffi- 
cient screening. The approximate T 2 behavior can, how- 
ever, also be obtained from the electron-phonon mech- 
anism with a reasonable distributions of coupling con- 
stants. Therefore we focus on this mechanism below. 
More recent experimental work (Vareka and Zettl, 1994) 
has suggested that p(T) at constant volume is linear in 
T down to relatively low T (100-200 K). This raises new 
questions, and the resistivity can not be considered as 
fully understood. 

The resistivity has been studied theoretically by sev- 
eral groups (Crespi et al, 1992; Erwin and Pickett, 1992; 
Antropov et al, 1993b; Lannin and Mitch, 1994), using 
the theory of Ziman (Grimvall, 1981). This theory pro- 
vides an approximate solution to the Boltzmann equation 
for the case of the electron-phonon scattering mechanism. 



P{T) = 



8tt 2 



f 

Jo 



u 2 p k B T 



cLu 



cosh(huj/k B T) - 1 



(10) 



where u. p is the plasma frequency. In the following it is 
assumed that the transport coupling function a^ r F(u) 
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can be replaced by the electron phonon coupling appro- 
priate for superconductivity 

a 2 F(w) = ^5^w„A„J(w-w„). (11) 

Due to the factor cosh(ftw/fcsT) — 1 in the denominator of 
Eq. (10), it is clear that phonons with a frequency much 
larger than T do not contribute, while they give a linear 
contribution in T if their frequency is much smaller than 
T. Using this approach, Crespi et al. (1992) demon- 
strated that theoretical results for the coupling constants 
to the H g modes alone of the type obtained by Varma et 
al. (1991) or Schlutcr et al. (1992) gave a poor fit of 
experiment, giving a much too weak temperature depen- 
dence at small T. The reason is that the H g phonons are 
too high in energy to give much contribution to p{T) at 
low T, in particular since in these theories most of the 
coupling occurs to the higher phonons. The couplings of 
Jishi and Dresselhaus (1992), emphasizing the lower H g 
phonons, gave better agreement. Crespi et al. (1992) 
argued that one also has to consider low-lying external 
modes, e.g., the intermolecular translational modes. In- 
cluding such modes, they showed that the couplings of 
Varma et al. (1991), Schluter et al. (1992) and Jishi 
and Dresselhaus (1992) give reasonable descriptions of 
the experimental T-dependence, if the coupling to the 
translational mode is chosen appropriately. 

As an example of such results, we show a calculation 
using the coupling constants deduced from photoemis- 
sion (Gunnarsson, 1995). In addition, the calculation 
included the coupling to the librations, using a theoret- 
ical coupling constant (Antropov et al, 1993b), and to 
the intermolecular modes, treating the coupling as an 
adjustable parameter, since its value is unknown. As in 
the calculation of Crespi et al. (1992), the energy of 
the intermolecular mode was set to 110 cm -1 , since neu- 
tron scattering shows many modes at this energy. Wc 
have set uj p = 1.36 eV. It was further assumed that some 
other mechanism, e.g., orientational disorder, provides a 
T-independent resistivity which accounts for the resis- 
tivity just above T c . The results are shown in Fig. 9. 
We can see that there is a good agreement with exper- 
iment over the whole temperature range. In particular, 
the approximate T 2 behavior is reproduced. This results 
from the contribution of many phonons at different ener- 
gies, whose contributions to the resistivity are gradually 
switched on one after the other as T is increased. For 
high temperatures in Fig. 9, the resistivity corresponds 
to a mean free path of the order of or less than the sepa- 
rations between the molecules. This raises serious ques- 
tions about the validity of the Boltzmann equation for 
these temperatures. 

The theoretical and experimental results in Fig. 9 are 
for constant pressure. The experimental data have also 
been converted to constant volume, by measuring the 



pressure dependence (Vareka and Zettl, 1994). The re- 
sistivity then shows an approximately linear behavior for 
temperatures larger than about 100 to 200 K. This lin- 
ear behavior appears not to have been explained theoret- 
ically, and it is an interesting problem. 

It is of particular interest to consider the absolute value 
of the resistivity increase with T, e.g., Ap = rho{T = 
260) — p(T = 0), since that may give information about 
the magnitude of A. Many estimates of the resistivity, 
however, are made for low temperatures, thereby giving 
the residual resistivity and not Arho. Since the shape 
of p(T), but not its absolute magnitude, appears to be 
well reproducible (Xiang et al., 1992), we use the exper- 
imental ratio Ap/p(T = 0) two relate etimates of Arho 
and p(T = 0). Here p(T = 0) refers to the resistivity 
just above the superconductivity transition temperature. 
Thus we convert our theoretical estimate of Arho to a 
theoretical estimate of p(T = 0) = 0.47 mil cm. It is 
interesting that the resistivity from orientational disor- 
der, 0.32 mflcm (Gelfand and Lu, 1992b; 1993), is not 
much smaller. Direct measurements for thin films and 
single crystals gave p(T = 0) ~ 1.2 nrOcm (Palstra et 
al, 1994) and p(T = 0) = 0.5 mUcm (Hou et al, 1993). 
Since the direct measurements involve various uncertain- 
ties (possibly somewhat granular structure for the thin 
films and geometrically very small samples for the single 
crystals), it is also interesting to compare with indirect 
measurements of the resistivity. Xiang et al. (1993) 
deduced p(T = 0) = 0.12 mficm from the paraconduc- 
tivity (additional conductivity due to superconducting 
fluctuations just above T c ) and Crespi et al. (1992) de- 
duced p(T = 0) = 0.18 mflcm from the upper critical 
field. Degiorgi et al. (1992) obtained p(T = 0) = 0.77 
nrOcm from a Kramer-Kronig analysis of the reflectiv- 
ity in the normal state and Rotter et al. (1992) found 
p(T = 0) = 0.4 nrOcm from the ratio between the reflec- 
tivity in the normal and superconducting state. These 
results are summarized in Table IX. 

Comparison between the theoretical and experimental 
estimates shows that the theoretical estimate falls in the 
correct range. Given the large range of experimental es- 
timates it is, however, hard to assess the accuracy of the 
theoretical calculation. We also observe that the theoret- 
ical result was obtained from a set of coupling parameters 
which give very strong couplings to the low-lying modes, 
which leads to a large resistivity, while several other esti- 
mates of the couplings would give a substantially smaller 
resistivity. 

The residual resistivity is smaller than the Mott limit, 
which, due to the small conduction electron density, is 
rather large for these systems, and has been estimated to 
be 1.4-2.0 milcm (Palstra et al, 1994). For high temper- 
atures, the resistivity corresponds to a mean free path 
shorter than the separation between the molecules. In 
such a situation the resistivity is usually expected to sat- 
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urate at a maximum value (Fisk and Webb, 1976). Such a 
saturation is not observed for the doped Cgo compounds 
up to rather high temperatures (T < 500 K) (Hebard 
et al, 1993), but for still larger temperatures it is seen 
for Rb 3 C6o (Hou et al, 1995). The saturation value was 
estimated to p sat = 6 ± 3 mfi-cm with the correspond- 
ing mean free path l sat = 1 ± 0.5 A, significantly smaller 
than the separation between the molecules. The reason 
for this small mean free path appears not to be under- 
stood, and it is an interesting question how one should 
conceptionally understand such a short mean free path. 

TABLE IX. The resistivity p(0) in mficm just above the su- 
perconductivity transition according to different theoretical 
and experimental estimates. For details see text. 



p(0) Method 

0.47 Theory. p(260) - p(0) -> p(0) 

0.32 a Theory. Disorder 

1.2 b Direct. Film 

0.5 C Direct. Single crystal 

0.12 d Paraconductivity 

0.18° Upper critical field 

0.77 f Reflectivity. K-K 

0.4 s Reflectivity. S-N 



a Gelfand and Lu (1992b, 1993) 
b Palstra et al. (1994) 
c Hou et al. (1993) 
d Xiang et al. (1993) 
e Crespi et al. (1992) 
f Degiorgi et al. (1992) 
s Rotter et al. (1992) 



6. Comparison of electron-phonon coupling estimates 

We are now in the position of comparing different the- 
oretical (Table IV) and experimental (Table VIII) esti- 
mates of the coupling constants. We first of all observe 
that the quantity A/V(0) is obtained in the calculations 
and is deduced from photoemission, while Raman and 
neutron scattering gives an estimate of AiV(0). Our as- 
sumptions about V(0) (7.2 states/eV per spin for K 3 C6o) 
is therefore crucial for these comparisons. 

There are substantial differences between the different 
calculations of A, between the different experimental es- 
timates, and between theoretical and experimental val- 
ues. Experimentally, a strong coupling to some of the 
low-lying modes is deduced. For instance, the total cou- 
pling to the lowest two modes is experimentally several 
times larger than theoretically. Thus the experimental 
estimates are 0.037, 0.048 and 0.059 eV, according to 
neutron, Raman and photoemission experiments, respec- 
tively, while the theoretical estimates are 0.004, 0.015, 
0.009 and 0.011 eV. For the high-lying phonons, the cal- 
culations give a rather strong coupling in agreement with 



the results from photoemission but in contradiction with 
the results from Raman. It is not clear why there is such 
a large discrepancy between the photoemission and Ra- 
man results for these modes. 

If a larger value for V(0) is assumed, the discrepancy 
between the theoretical values and the Raman and neu- 
tron estimates for the low-lying modes would be reduced, 
but the discrepancy with the photoemission values would 
remain. At the same time the difference between the 
theoretical and the Raman estimates would be increased 
for the high-lying modes. As discussed above, relatively 
small errors in the phonon eigenvectors could lead to a 
large underestimate of the theoretical coupling to the 
low-energy modes. 



B. Alkali phonons 

Very early it was pointed out that the alkali modes 
may provide a strong electron-phonon coupling (Zhang 
et al, 1991a). Each Ceo molecule in A 3 Ceo is sur- 
rounded by 14 alkali atoms (6 octahedral and 8 tetra- 
hedral). The corresponding phonon frequencies are low 
(0.013 eV (Prassides et al, 1992); 0.014 eV (Renker et 
al, 1993); 0.13 eV (Mitch and Lannin, 1993)), implying 
weak force constants. If an electron is transferred to a 
Ceo molecule, one may then expect the surrounding alkali 
ions to move towards the Cgo molecule, giving an efficient 
screening. It has been estimated that the alkali phonons 
could contribute an attractive interaction of the order 0.9 
eV (Zhang et al., f 99f a), which corresponds to A ~ 2, i.e., 
a strong coupling. Later experiments demonstrated that 
T c is approximately independent of the alkali atom mass 
(see Sec. IIA and IIC), showing that the alkali modes 
must contribute very little to T c . This was rationalized 
in terms of an efficient metallic screening (Gunnarsson 
and Zwicknagl, 1992). If an electron is transferred to a 
Ceo molecule, there is a tendency to screen it by transfer- 
ring charge from that molecule to more distant molecules. 
The alkali ions then couple essentially to the net charge 
change of the Cgo molecule. Within the RPA this screen- 
ing mechanism is very efficient, and when an electron is 
transferred to a Cgo molecule the change of its net charge 
is only a few hundredths of an electron (Gunnarsson and 
Zwicknagl, f992). The coupling to the alkali modes is 
then reduced by one to two orders of magnitude. Al- 
though RPA may greatly overestimate the efficiency of 
the screening, this nevertheless suggests a large reduc- 
tion of the coupling. 

The alkali phonon frequency for the tetrahedral po- 
sition has been calculated using the density functional 
formalism in the LD approximation by Bohncn et al. 
(1995). They found the energy 0.015 eV, in good agree- 
ment with experiment. 
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C. Librations 

Early point contact measurements of the superconduc- 
tivity gap gave a large value of 2A/T C (Zhang et al., 
1991b), suggesting strong coupling effects and the cou- 
pling to low energy modes. It was therefore suggested 
that the coupling to the librations may be of importance 
(Dolgov and Mazin, 1992; Mazin et al, 1993b). There 
are several factors influencing this coupling. For the in- 
tramolecular phonons, the main changes in the matrix 
elements come from changes in bond lengths and bond 
angels, i.e., changes of the order of the intramolecular 
energy scale. For the librations, which involve rigid ro- 
tations of the molecules, only the hopping between the 
molecules is influenced, i.e., changes of the order of the 
intermolccular energy. The amplitude factor in Eq. (6) 
also tends to be small for the librations due to the large 
moment of inertia, la ~ 4QMR 2 , where M is the carbon 
mass and R is the Cqq radius. On the other hand, the 
very small libration frequency (~ 4 meV, (Christides et 
al, 1992 )) tends to increase both the amplitude factor 
in Eq. (6) and the A in Eq. (5). It is therefore a priori 
not easy to estimate the magnitude of the coupling to 
the librations. 

We observe that the simplifications in the calculation 
of the coupling to the intramolecular phonons (Lannoo et 
al, 1991) do not apply to the librations and the Brillouin 
zone sums in Eq. (5) should be performed. This has been 
performed within a tight-binding formalism, which gave 
Xi/N(0) ~ 0.01 (Antropov et al., 1993b). Comparing 
the splitting of the t lu band at the T point in the tight- 
binding calculation and a ASA-LMTO LDA calculation, 
suggested that A; should be reduced by about a factor 
10, giving Xi/N{0) ~ 0.001 and A; - 0.01 (Antropov et 
al, 1993b). Thus the coupling to the librations is found 
to be very weak. 

The coupling has also been studied experimentally 
(Christides et al, 1992; Rcznik et al, 1994). They stud- 
ied the librations above and below T c . Since the libration 
energy is smaller than the superconductivity gap 2A, the 
broadening of the libration due to the decay in electron- 
hole pairs should go to zero for small temperatures (Zey- 
her and Zwicknagl, 1990). Experimentally, however, no 
change in the width was observed when the temperature 
was lowered below T c . It was therefore concluded that 
A; < 0.08, since the change in width should otherwise 
have been observed (Christides et al., 1992). 

D. Intermolecular modes 

The lattice-dynamics of K 3 Ceo has been studied by Be- 
losludov and Shpakov (1992). They included the van der 
Waals interaction between the Cgo molecules, the elec- 
trostatic interactions between the ions, Born-Mayer type 
of repulsive interactions to describe the repulsion of the 



ion cores and a valence force model to describe the cova- 
lent interaction between the carbon atoms inside a Cgo 
molecule. The parameters were adjusted to reproduce 
the experimental intramolecular T\ u vibrations. This 
lead to acoustic modes which at larger wave vectors are 
mixed with librational modes at about 3.5 meV. At larger 
energies (6 meV and 16 meV) they found bands of al- 
kali and Ceo character. This seems to be consistent with 
neutron scattering experiments (Prassides et al, 1991; 
Christides et al., 1992). The translational phonons have 
also been studied by Zhang et al. (1992) for K 3 Ceo and 
by Wang et al. (1991) for C 60 . 

Pickett et al. (1994) have studied the electron-phonon 
coupling for translational modes. They used a tight- 
binding parameterization of the hopping integrals. A 
distortion of the lattice corresponding to the longitudinal 
zone boundary (110) intermolecular optic mode was in- 
troduced. From the displacement of the band structure 
energies they estimated that A ~ 0.01 for these modes. 
The intermolecular modes should therefore have a small 
effect on T c . 

IV. COULOMB INTERACTION AND METALLICITY 

In this section we discuss the Coulomb interaction be- 
tween two electrons on a free Ceo molecule ((To) and for a 
molecule in the solid (U). The large Coulomb interaction 
has played an important role in the discussion about su- 
perconductivity in the fulleridcs. This has ranged from 
the issue of why A3C60 are not Mott-Hubbard insula- 
tors (Lof et al, 1992), over questions about whether or 
not retardation effects can reduce the Coulomb inter- 
action enough to allow for superconductivity driven by 
the electron-phonon interaction (Anderson, 1991; Gun- 
narsson and Zwicknagl, 1992) to the suggestion that the 
Coulomb interaction could actually be the cause of su- 
perconductivity (Chakravarty et al, 1991; Baskaran and 
Tosatti, 1991; Friedberg et al, 1992). 

A. Coulomb Uq for a free molecule 

A simple estimate of Uq for a free molecule can be ob- 
tained by assuming that the charge from an electron is 
uniformly spread out over a sheet with the radius R. The 
interaction between two electrons is then e 2 /R = 4.1 eV, 
if R = 3.5 A. This estimate neglects the relaxation of the 
charge density when an electron is added to the molecule. 
To include this, one has to perform self-consistent calcu- 
lations for different number of electrons. This has been 
performed by several groups. Most of these calculations 
were based on the local density approximation of the den- 
sity functional formalism (Hohcnbcrg and Kohn, 1964; 
Kohn and Sham, 1965; Jones and Gunnarsson, 1989). 
Uq is given by the change of the position of the t\ u level 
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when charge is added to this level. The calculated values 
were 3.0 (Pederson and Quong, 1992), 2.7 (Antropov et 
al, 1992), 3.0 (de Coulon et al, 1992) and 3.1 eV (Martin 
and Ritchie, 1993). 

Estimates of Uq can also be obtained from experimen- 
tal data. One estimate is provided by 

U = /(C 60 ) - A(C 60 ) -E g = 3.3 eV, (12) 

where I(Cqo) — 7.6 eV (de Vries et al, 1992) is the ion- 
ization energy, A(C^) = 2.7 eV (Hcttich et al, 1991) is 
the affinity energy and E g — 1.6 eV (Gensterblum et al, 
1991) is the band gap as measured by photoabsorption 
for a free Cgo molecule. An alternative measure is given 

by 

= /(Ceo) - MC m ) ~ 2.7 eV, (13) 

where I(Cq ) = 2.7 eV (Hettich et al, 1991) is the ion- 
ization energy and and A(Cq ) is the affinity energy of 
CgQ. Cqq has been observed experimentally (Hcttich et 
al., 1991; Limbach et al., 1991) and if it is indeed sta- 
ble, A(Cqq) > and Uq < 2.7 eV. It may however, be 
metastable with A(Cg ) slightly negative and Uq slightly 
larger than 2.7 eV. The estimate (Eq. (12)) refers to the 
interaction between a hole in the h u level and an elec- 
tron in the t\ u level, while the estimate (Eq. (13)) refers 
to the interaction between two electrons added to a Ceo 
molecule. Since the electron and hole try to stay together 
while the two electrons try to stay apart, it is not sur- 
prising that the latter estimate is smaller (Antropov et 
al, 1992). 



B. Coulomb U for the solid 

In the solid, there are additional screening effects due 
to the polarization of the surrounding molecules when 
two electrons are added to a given molecule. In the fol- 
lowing, we have a Hubbard model in mind, where the 
on-site Coulomb interaction is included, but the screen- 
ing due to the polarization of the surrounding molecules 
is only taken into account implicitely as a renormalization 
of an effective U. In this way the Coulomb interaction 
for a solid was found to be U = 1.27 eV (Pederson and 
Quong, 1992) and 0.8-1.3 eV (Antropov et al, 1992). The 
Coulomb interaction may for many purposes be drasti- 
cally reduced by the metallic screening. This reduction 
should, however, not be included if U is used in, say, a 
Hubbard model, since the metallic screening is then in- 
cluded explicitly in the model. To reduce the parameter 
U by the metallic screening would then involve double 
counting. 

Experimentally, U has been estimated from Auger 
spectroscopy. Without any Coulomb interaction the 



Auger spectrum would be just a convolution of the one- 
particle spectrum, while the presence of an on-site molec- 
ular interaction U should displace the spectrum by ap- 
proximately U. The comparison of the Auger spectrum 
with the convoluted photoemission spectrum then gives 
an estimate of U = 1.6 eV (Lof et al, 1992) and U = 1.4 
eV (Briihwiler et al, 1992). For the highest occupied or- 
bital (h u ) the estimate of U is about 0.2 eV smaller (Lof 
et al, 1992), and for the t\ u orbital it may be still some- 
what smaller. Due to surface effects, Auger spectroscopy 
may also slightly overestimate U (Antropov et al, 1992). 

C. Mott-Hubbard transition in systems with a large 
degeneracy 

From the estimates of U and W it follows that the ratio 
U /W is substantially larger than one. It is then gener- 
ally expected that the system should be a Mott-Hubbard 
insulator (Georges et al, 1996). It has therefore been 
argued that stoichiometric A 3 Ceo must be an insulator 
and that the experimental samples of A 3 Ceo are metal- 
lic only because they are nonstoichiometric (Lof et al, 
1992). This would be a similar situation as for the high 
T c superconductors, and it might suggest that various 
properties of A3C60 should be understood in a similar 
way as for the high T c supersoncutors. We therefore dis- 
cuss the origin of the metallicity in A3C60 in more detail. 

The estimates of the critical value of U/W for a Mott- 
Hubbard transition have been made mostly for systems 
without orbital degeneracy, while the t\ u orbital has a 
threefold degeneracy. This was considered in early work 
of Lu (1994), who used the socalled Gutzwiller Ansatz 
and approximation, finding that the orbital degeneracy 
leads to a Mott-Hubbard transition for U/W <~ (N + lj, 
where N is the orbital degeneracy. This strong degen- 
eracy dependence raised questions about its qualitative 
explanation and about why not almost all systems with 
partly filled bands are metals. 

To discuss qualitatively the effect of orbital degeneracy, 
we consider a half-filled system with the orbital degener- 
acy N and M sites. The band gap is then given by 

E g = E(K+l) + E(K-l)-2E(K), (14) 

where K = NM is the number of electrons in the neu- 
tral state and E(K) is the ground-state energy for K 
electrons. If we consider the system with K electrons in 
the limit of a large U, configurations with N electrons 
per site dominate the wave function. The hopping of an 
electron costs an energy U, and it is therefore strongly 
suppressed. In the state with K + 1 electrons, the ex- 
tra electron can, however, hop without an extra cost in 
Coulomb energy The extra occupancy can move through 
the hop of any of N electrons, which gives an extra fac- 
tor \/N in the corresponding hopping matrix element 
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compared with the one-electron case (Gunnarsson et al., 
1996). This suggests that the energy is lowered by about 
VNeb, where et, is the bottom of the one-particle band. 
Similar arguments for the K — 1-electron state, suggest 
that 



E„ 



U 



NW. 



(15) 



Exact diagonalization calculations and Monte Carlo cal- 
culations for A3C60 (Gunnarsson et al., 1996) support 
these conclusions and suggest that they can be extrapo- 
lated to intermediate values of U . Detailed calculations 
lead to the conclusion that the Mott-Hubbard transition 
takes place for U/W ^ 2^. This is at the upper limit 
of both experimental and theoretical estimates of U /W. 
This suggests that A3C60 is on the metallic side of a 
Mott-Hubbard transition. Although experimental sam- 
ples of A3C60 may be nonstoichiometric, there is then 
no need to assume so to explain why they are metallic. 
Since almost all physical systems have an orbital degen- 
eracy, the degeneracy dependence of the Mott-Hubbard 
transition should be of interest for many other systems. 



A. Retardation effects. Summing ladder diagrams 



It was very early pointed out that the Coulomb re- 
pulsion can be reduced by retardation effects (Bogoli- 
ubov et al., 1958; Morel and Anderson, 1962; Schrieffer, 
1964; Ginzburg and Kirzhnits, 1982). Normally, the elec- 
tronic energy scale is much larger than the energy scale 
of the phonons. By undergoing multiple scattering into 
higher energy states, the electrons can move in a corre- 
lated way and thereby reduce the Coulomb interaction 
but still have an attractive interaction via the phonons. 
This is shown schematically in Fig. 10 for the scattering 
up into a higher subband. Mathematically, this is usually 
expressed by summing ladder diagrams (see Fig. 11). In 
simple models, e.g., assuming that all Coulomb matrix 
elements are equal, the ladder diagrams can be summed, 
giving the effective interaction (Bogoliubov et al., 1958; 
Morel and Anderson, 1962; Schrieffer, 1964; Ginzburg 
and Kirzhnits, 1982) 



V. COULOMB PSEUDOPOTENTIAL 

The electron-phonon interaction provides an attrac- 
tive interaction between the electrons, which may be of 
the order 1/10 eV. This interaction is counteracted by 
the Coulomb repulsion, and two electrons on the same 
molecule in an undoped Cgo solid may feel a repulsion of 
the order U ~ 1 — | eV (see Sec. IVB). In this section we 
first discuss how the Coulomb repulsion may be reduced 
due to screening and retardation effects. A dimensionless 
quantity 



li = UN(0) 



(16) 



is introduced, where U is some typical (screened) inter- 
action and iV(0) is the density of states per spin. Due to 
retardation and other effects, [i is renormalized to \i* , the 
Coulomb pseudopotential. A closely related issue is the 
question if the Coulomb interaction could be so effec- 
tively reduced that it becomes attractive, providing an 
electronic mechanism for superconductivity. This goes 
back to calculations for a model of a free molecule by 
Chakravarty et al. (1991), and is discussed in the last 
part of this section. 



1 + fAog(B/uj ph ) ' 



(17) 



where B is a typical electron energy (half the band width) 
and LUph is a typical phonon energy. If B/u p h >> 
1, /i* can be strongly reduced relative to yu. In the 
limit fAn(B/u; ph ) » 1, Eq. (17) simplifies to ^* « 
l/log(B/uj p f l ), which may be of the order 0.1-0.2. The 
importance of this reduction is illustrated by the McMil- 
lan formula (McMillan, 1968) for the transition temper- 
ature 



T ^iexpf L ° 4 ( 1 + A ) 1 

c ~ 1.2 P[ A-/x*(l + 0.62A) J 



In this formula, the relevant quantity is 



(18) 



1 + 0.62A A 



(1 + A) 



1 + A (1 + A) 



- H 



(19) 



It is immediately clear that [i* is as important as A for 
the T c . 



FIG. 10. Schematic picture of two subbands and the virtual 
Coulomb scattering of two electrons into the higher subband. 



> 1 > 



> i > 1 > 
> 1 > 1 > 



FIG. 11. Ladder diagrams describing the repeated scatter- 
ing of two electrons. A full line represents an electron and a 
dashed line the (screened) Coulomb interaction. 
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Doped Ceo compounds have a number of narrow sub- 
bands with a width of about 1/2 eV spread out over 
about 30 eV (see Fig. 2). It has been argued that the 
relevant energy scale is the total width of all the sub- 
bands (Varma, 1991). Using Eq. (17) one then finds that 
li* <~ 0.2 is only slightly larger than what is believed to 
be the case for conventional superconductors. Together 
with current estimates of A, values of T c of the right order 
of magnitude can then easily be obtained (Varma et al, 
1991; Schluter et al, 1992; Mazin et al, 1992). On the 
other hand, Anderson (1991) asserted that the relevant 
energy scale is the width of the t iu subband. Since this 
energy is comparable to the phonon energies, the retar- 
dation effects are then expected to be small. Anderson 
therefore argued that [i* is large, and that the phonon 
mechanism alone cannot be enough for explaining the 
superconductivity. 

We first discuss the retardation effects on n* within 
the formalism described above, where ladder diagrams in 
the screened Coulomb interaction are summed. Later we 
discuss effects beyond this formalism. A crucial quan- 
tity is the Coulomb matrix element for scattering two 
electrons from the t\ u band into higher subbands (see 
Fig. 10). If these matrix elements are very small, the 
higher subbands should play a small role in this for- 
malism. Chakravarty et al. (1992b) analyzed this in 
a model where they included the on-site interaction Uc 
between two electrons on the same carbon atom and the 
interaction Vc between two electrons on different carbon 
atoms on the same Cgo molecule. They then found that 
the Coulomb scattering matrix element between differ- 
ent subbands is of the order Uc/Q0, while the interaction 
within the same band is Uc/Q0 + Vc- Assuming that 
U c « 5 - 10 eV and V c « 0.5 - 1.0 eV, they found that 
the interband interaction is very much smaller than the 
intraband interaction. 

The matrix elements of the screened Coulomb inter- 
action have been studied in the RPA (Gunnarsson and 
Zwicknagl, 1992; Gunnarsson et al., 1992). The sys- 
tem was described in a tight-binding approach (Satpa- 
thy, 1992; Laouini et al, 1995) with one 2s- and three 
2p-orbitals per carbon atom, and the screening was cal- 
culated with local-field effects included. A long-range 
Coulomb interaction was included, which involved an in- 
teraction between different atoms both on the same Ceo 
molecule and on different molecules. The Coulomb in- 
teraction on a carbon atom was set to Uc = 12 eV. A 
screening calculation within the LDA approximation has 
been performed by D.P. Joubert (1993), with fairly simi- 
lar results. The tight-binding results are shown in Table 
X. The table shows both intraband and interband inter- 
actions 

TF(a/3, 7 <S,q) = (20) 
N [ d 3 rd 3 rV; k (r)^ k+q (r)VF(r,r> 7k (r )^ k+q (r ), 



where 



is a Bloch sum of molecular orbitals ip a (r), W(r,r ) is 
the screened interaction and TV is the number of unit 
cells. The intraband integrals correspond to the interac- 
tion of charges which integrate to one in the unit cell, 
i.e., these integrals are dominated by a monopole inter- 
action. On the other hand, (/) Q (r)0 / g(r) with a ^ j3 in- 
tegrates to zero, and therefore the interband matrix el- 
ements W (a/3, a/3) correspond to multipole interactions. 
In the solid and without screening, the intraband ma- 
trix elements are therefore very large for small q, since 
they represent the interaction of monopoles placed on 
the different molecules and with slowly varying phase 
factors. When screening is introduced in the undoped 
solid, these integrals are reduced due to the polariza- 
tion of the molecules. This reduction is much more dra- 
matic in the doped solid due to the effective metallic 
screening, which involves a charge transfer between the 
molecules. The unscreened interband matrix elements 
are much smaller than the unscreened intraband ele- 
ments, due to the shorter range of the multipole interac- 
tion. In the undoped solid they are further reduced, pri- 
marily due to rearrangement of the charge density inside 
the molecule. In the doped solid, the metallic screening 
does not essentially further reduce these matrix elements, 
since a charge transfer between the molecules does not 
contribute to the screening of multipoles. The moderate 
reduction shown in Table X is due to the change in the 
filling of the levels and the new channels (in particular, 
tiu — > t\g) opened up for intramolecular screening. It 
has been argued that the intraband (metallic) screening 
should not be included in the calculation of the interband 
matrix elements (Lammert et al, 1995). The calcula- 
tions discussed above show, however, that the metallic 
screening of the interband matrix elements is negligible 
(Gunnarsson et al, 1992). 

The Coulomb integral W (a/3, a/3) with ct / |3 refer- 
ring to the t\ u orbitals is also of great interest, since it 
may tend to favor singlet states (Lammert and Rokshar, 
1993). With only the short-range on-site carbon inter- 
action Uc, it is Uc/W0 — 0.12 eV for a free molecule 
(Lammert and Rokshar, 1993). Inclusion of the long- 
range Coulomb interaction reduces the integral to 0.09 
eV, implying that the long range of the Coulomb inter- 
action is not very important for this integral (in con- 
trast to the monopole integral W(aa, act)). Allowing for 
screening W(a(3, a/3) is reduced to 0.023 eV. This value is 
essentially unchanged if the screening by the surrounding 
molecules is included. 
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TABLE X. The intraband Coulomb interactions (Eq. 
(20))between two equal W{aa,aa,c\) and two different 
W(aa, (3(3, q) ti u conduction band Bloch states as well as the 
largest interband scattering matrix element W tnter (ti u — > 
We consider both a free molecule (Mol) and a solid 
(Sol), where in the latter case q = G/2 and q = G/20, where 
G = 27r/o(l, 1, 1). The unscreened (Unscr) and screened 
(Scr) interactions in the undoped (Undop) and, for the solid, 
doped systems are shown. 



Svs 



'(aa, aa, q) 



W 



/?/?,q) 



W 



G/20 G/2 G/20 G/2 G/20 G/2 



Unscr Mol 

Scr Mol 

Unscr Sol 

Undop Sol 

Doped Sol 



3.81 

3.72 
176.8 2.03 
78.9 1.68 
0.094 0.065 



3.63 

3.67 
176.7 1.86 
78.8 1.63 
0.060 0.031 



0.77 

0.14 
1.01 0.98 
0.15 0.15 
0.10 0.10 



Because of the very efficient metallic screening of the 
intraband matrix elements, at least within RPA, these 
matrix elements are smaller than the largest interband 
matrix elements. These results depend crucially on the 
fact that these matrix elements have very different char- 
acter (monopole versus multipole interaction). This be- 
havior is also very different from the electron gas, since 
there is an abrupt change in character when going from 
one subband to another, even if they are close in energy, 
due to the fact that each subband is derived from a dif- 
ferent molecular level with a different spatial character. 
Since Eq. (17) is derived assuming that all matrix ele- 
ments are equal, it can not be used here. Approximate 
summing of the ladder diagrams, using the calculated 
matrix elements, shows that n* is close to zero in this 
formalism (Gunnarsson and Zwicknagl, 1992). 

It is interesting to ask why the intraband matrix el- 
ements are so small. The screened interaction can be 
written as 



w = (i-vpy 



(22) 



where W is a matrix representing the different Coulomb 
matrix elements, v is the corresponding unscreened ma- 
trix elements and P is the polarizability, which is here cal- 
culated in the RPA. For the intraband matrix elements, 
the tiu orbitals dominate the screening, and we can trans- 
form to the corresponding basis, neglecting all other or- 
bitals. For small |q|, the diagonal elements of P are then 
related to the density of states N(0) per spin. Since both 
N(0) and the diagonal elements of v are large, we can 
assume that the product is much larger than unity. Es- 
sentially, the two factors v in Eq. (22) then drop out, and 
with appropriate assumptions one can then derive (Gun- 
narsson et al, 1992) that the intraband matrix elements 
are 

^W(aa, aa, q) + ^W(aa, (3(3, q) w , (23) 



for small |q|. To obtain the unrenormalized /i we multi- 
ply by N(0) obtaining /i w 0.5. Averaging over q gives 
fi rts 0.4. We note that W becomes very small in the RPA 
for large N(0), since in the RPA the cost of screening is 
purely a kinetic energy cost, and for a large N(0) this 
cost is very small. For large values of U, we expect RPA 
to become a poor approximation, and it is therefore an 
interesting and important question how accurate RPA is 
in this case. 



B. Beyond ladder diagrams 

The considerations above were all within a framework 
where ladder diagrams in the screened interaction arc 
summed. This raises the question about the validity 
of neglecting other diagrams. For instance, Grabowski 
and Sham (1984) studied the lowest order vertex correc- 
tions for the electron gas and found important correc- 
tions. Gunnarsson et al. (1992) constructed other sets 
of diagrams, which to a large extent cancelled the ladder 
diagrams. This raises important questions about which 
diagrams to include. 

To circumvent this problem, Gunnarsson and Zwick- 
nagl (1992) and Gunnarsson et al. (1992) studied a two- 
band model where some exact results can be obtained. 
Thus they considered the Hamiltonian 

2 2 

H = X XX £nUm ° +t EEE ^Lvw 

i n—1 a <\?> n—1 a 



1 nn aa 

+^12X^2^121^11^11 + h.c], 



(24) 



where i is a site index, e n the energies of the two levels, 
t is the hopping integral, Uu is the intraband interac- 
tion and U12 is the interband scattering. The particu- 
larly strong interband scattering between the t\ u and t\ g 
bands in Ceo is of a similar type, namely the excitation 
of two electrons between two bands (see Fig. 10). It was 
then assumed that 



Ae = e 2 - £1 >> U\2,t. 



(25) 



In this limit, the upper subband can be projected out 
and a new effective Hamiltonian be obtained, which de- 
scribes low-energy properties correctly. This process does 
not generate any new terms for the Hamiltonian (24), but 
it renormalizcs the intraband interaction to 



Ueff = U U 



2Ae' 



(26) 



If it is now assumed that U~i 2 /(2Ae) « Un, the renor- 
malization of the intraband interaction is small. The 
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properties of this effective Hamiltonian then differ very 
little from the one-band model, where the upper band 
was completely neglected. Thus the upper band has a 
small influence on the low energy properties in general 
and n* in particular. For instance, in RPA the screened 
interaction is 1/[27V(0)] in both the one-band and two- 
band models, if U X1 N{Q) » 1. 

On the other hand, the ladder diagrams in the screened 
interaction due to scattering into the upper subband can 
be summed. If it is assumed that 



n 2 
< 12 



N(0) ~ 2Ae 



(27) 



this approach predicts a large renormalization of fi* . The 
ladder diagrams subtract a quantity C/f 2 /(2Ae), which, 
although small compared with U\\, is large compared 
with the screened interaction in the limit considered. 
Thus the summation of ladder diagrams leads to a qual- 
itatively incorrect result in this case. 

The projection method and the summation of ladder 
diagrams describe similar physics. The difference is that 
in the rigorous projection method, the high energy de- 
grees of freedom were eliminated first, while in the lad- 
der diagram approach the low energy degrees of freedom 
were treated first (by introducing the metallic screening) 
and the high energy degrees of freedom later (by sum- 
ming the ladder diagrams). These results should be of 
relevance for the retardation effects due to the high-lying 
states also for conventional superconductors. 

The model above contains important features of A3C60 
and the assumptions are not too unrealistic for this sys- 
tem. The results therefore suggest that the summation of 
ladder diagrams in the screened interaction gives unreli- 
able results for A3C60 and that the renormalization of fj,* 
due to the higher subbands may not be large. For simplic- 
ity one may therefore assume that this renormalization 
can be neglected completely. If in addition RPA screen- 
ing is used, one arrives at ji* <~ 0.4. The approach above 
has only addressed the renormalization due to the higher 
subbands but not due to processes within the t\ u band 
itself. If we assume that the traditional theory is valid 
for these processes, /i* is renormalized to 0.3 according 
to formula (17), if we set B = 0.25 eV and u p h = 0.1 eV. 
Such a value is then to be used in the McMillan formula, 
while the value 0.4 applies to the Eliashbcrg equation in- 
cluding one subband. We notice that only static screen- 
ing was considered above, although inclusion of dynamic 
effects may influence the results substantially (Rietschel 
and Sham, 1983). 

A very interesting work was performed by Takada 
(1993). He considered a model including electrons in 
the ti u band, phonons and plasmons, where the plas- 
mons screened the strong Coulomb repulsion. Guided by 
the Ward (1950) identity and work by Nambu (1960), 
Takada used an approximate vertex correction. In this 



way he was able to include corrections to Migdal's the- 
orem, vertex corrections to the electron-electron scatter- 
ing in the t\ u band and the dynamical screening of the 
Coulomb interaction. It was found that the corrections 
to Migdal's theorem enhanced the value of T c for moder- 
ate values of the phonon frequency. Furthermore it was 
found that the vertex corrections to the electron-electron 
scattering reduce T c . Thus there is a partial cancelation 
of the enhancement (Rietschel and Sham, 1983) of T c 
due to dynamic screening effects, as had also been found 
earlier (Grabowski and Sham, 1984). Takada used the 
value [i = 0.4 for the Coulomb pseudopotential derived 
above, to describe static screening effects and the renor- 
malization due to scattering in other sub bands than the 
ti u band. Using A = 0.58 and adjusting the electron- 
plasmon coupling, Takada (1993) could then reproduce 
the experimental T c for K 3 C6o- The corresponding iso- 
tope effect a = 0.15 for K3C60 was, however, surprisingly 
small. To describe the lattice parameter dependence of 
T c , Takada had to assume that [i = 0AW o /W, where W 
and W are the band widths of K3C60 and the system 
studied, respectively. This is somewhat surprising, since 
in RPA [i is predicted to be essentially independent of 
the lattice parameter. 



C. Electronic mechanisms 

Chakravarty and Kivelson (1991), Chakravarty et al. 
(1991) and Baskaran and Tosatti (1991) proposed that 
there may be an effective attractive interaction between 
the electrons for a Ceo molecule of a purely electronic ori- 
gin. Chakravarty and Kivelson (1991) introduced a Hub- 
bard model of the 7r-electron system in a free molecule 
with an on-site interaction Uc, only. Using second order 
perturbation theory, they calculated the energy of Ceo, 
C^q and CgQ~ . This provided an estimate of the effective 
repulsion for CgQ. When the singlet state of C^q was con- 
sidered, the result consisted of a positive, linear term and 
a negative, quadratic term in Uc- The interaction there- 
fore becomes attractive for sufficiently large Uc ~ 3io, 
where to is the intramolecular hopping. The tendency to 
the formation of a singlet state results from an integral 
of the type W(a/3,a/3) (see Eq. (20)) with a ^ (3 be- 
longing to the ti u band (Lammert and Rokhsar, 1993), 
which may become negative when the higher states are 
projected out in second order perturbation theory. The 
following discussion has focusscd on the validity of the 
second order perturbation theory and the neglect of a 
long-range Coulomb interaction in this type of work. 

White et al. (1992) addressed the question about the 
validity of second order perturbation theory, by studying 
small model systems. In particular, they considered a 
model with twelve atoms located on a truncated tetra- 
hedron. For this system exact diagonalization can be 
performed and compared with second order perturbation 
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theory. White et al. (1992) concluded that second order 
perturbation theory is qualitatively correct for interme- 
diate values of Uc in the sense that exact diagonaliza- 
tion also gives an attractive interaction, but that second 
order perturbation theory may greatly overestimate the 
magnitude. Auerbach and Murthy (1992) and Murthy 
and Auerbach (1992) studied a model with electrons on 
a sphere and with an interaction of variable range. They 
calculated the second order contribution to pairing and 
they also estimated the third order terms. The second 
and third order contributions were found to be compa- 
rable at an interaction strength which was substantially 
smaller than the strength giving pairing in the second 
order calculation. This again raises doubts about the 
accuracy of the second order calculation. Berdenis and 
Murthy (1995) performed a rcnormalization group cal- 
culation for a Hubbard model of Ceo using a method of 
Tokuyasu et al. (1993). This method remains accurate 
up to larger values of the atomic Uc than the second or- 
der perturbation theory, but its accuracy becomes ques- 
tionable at roughly the value of Uc where it predicts an 
attractive interaction. Thus it seems to still be unclear 
whether or not the interaction can be attractive for an 
on-site Hubbard model of a free Ceo molecule. 

Goff and Phillips (1992; 1993) and Auerbach and 
Murthy (1992) have studied the effective interaction as 
a function of the range of the interaction. They found 
that the use of a more long-ranged interaction tended to 
suppress the formation of an attractive effective interac- 
tion. With an on-site interaction, two electrons added 
to the t\ u orbitals can almost avoid being on the same 
site simultaneously and they have only a weak direct in- 
teraction W(aa,aa). It is then not entirely surprising 
that the renormalized interaction W (a/3, a/3) can lead to 
an effective attractive interaction. If the interaction be- 
comes more long-ranged, the direct interaction rapidly 
increases, leading to a repulsive effective interaction. For 
the free molecule, the interaction is long-range and there 
is actually an "anti-screening" effect for electrons on the 
opposite side of the molecule (Gunnarsson et al, 1992). 
The reason is that for a finite system, screening charge 
cannot be moved to infinity as for an infinite system. It 
is therefore no surprise that the effective interaction for 
the free molecule is large and repulsive, as discussed in 
Sec. IV. A. Chakravarty et al. (1991) therefore from the 
start pointed out that the attractive electronic mecha- 
nism could become operative only if the screening inside 
the molecule is much more efficient in the metal than in 
the free molecule. 

Lammert et al. (1995) considered a model with a 
long-range Coulomb interaction. It was argued that 
the molecules surrounding a specific molecule provide 
a metallic environment, which was approximated by a 
metal with a spherical cavity with the radius R c . In this 
way the corresponding degrees of freedom were integrated 



out first, although the metallic screening corresponds to 
a low energy scale in the problem (W ~ w p i <~ 0.5 eV 
(Knupfer et al, 1995)). The resulting screened inter- 
action was calculated and used for the low-frequency in- 
teractions in the problem. The effective singlet two-body 
interaction was then calculated diagrammatically, includ- 
ing the the first and second order diagrams. An attractive 
interaction was obtained for R c in the range 5-6 A and 
the on-site carbon Uc in the range 8-10 eV. 

An alternative way of treating the problem within this 
spirit, is to express the metallically screened interaction 
Umet.scr. in terms of the Coulomb interaction J7 between 
two electrons on a free Ceo molecule. It is first assumed 
that in the solid the addition of a charge q to a molecule 
induces a potential —qVo due to the metallic screening 
by the surrounding molecules, reducing [7 to 

Umet.scr. = Uq — Vq (28) 

in the solid. We here assume that due to the small 
hopping between the Ceo molecules, the surrounding 
molecules respond to a spherical symmetric average of the 
charge on the central molecule. In the model of Lammert 
et al. (1995), replacing the surrounding by a spherical 
cavity, Vo = e 2 /R c . To obtain an estimate of Vo (or R c ) 
in the above spirit, one can treat the surrounding Cqq 
molecules as grounded, metallic spheres with the radius 
4.46 A (which reproduces the experimental molecular po- 
larizability) . This leads to Vo = 2.2 eV (Krier and Gun- 
narsson, 1995), corresponding to R c = 6.6 A, i.e., too 
large to obtain an attractive interaction in the work of 
Lammert et al. (1995). If this is combined with the 
estimate of Uq <~ 2.7 eV (see Sec. IVA), a large and pos- 
itive U me t.scT. ~ 0.5 eV (fj, = N(0)U ~ 4) results. This 
may suggest that it is important to also allow for charge 
transfer to the molecule where the studied electrons are 
located, as normally assumed for metallic systems. 

An alternative electronic mechanism was suggested by 
Friedberg et al. (1992). They argued that it is favorable 
to form a two-electron wave-function involving one t\ u 
and one t\ g state. In this way a correlated state can be 
formed, which lowers the Coulomb energy. Friedberg et 
al. (1992) argued that this more than outweighs the cost 
in one-electron energy. This two-electron state is boson- 
like and forms a Bose band. It was suggested that this 
may lead to a Bose-Einstein condensate. 

VI. BEYOND THE MIGDAL-ELI ASH BERG THEORY 
A. Migdal's theorem 

Migdal's theorem states that vertex corrections in the 
electron-phonon interaction can be neglected if the typi- 
cal phonon frequencies u> p h are sufficiently much smaller 
than the electronic energy scale, say the Fermi energy 
Ep (Migdal, 1958). This theorem is assumed to be valid 
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in the derivation of the Eliashbcrg (1960) theory. For 
the fullerides, the highest phonon frequency is about 0.2 
cV and the width of the t\ u band is about 1/2 eV, giv- 
ing Ep ~ 0.25 eV. Migdal's theorem may therefore be 
strongly violated for the fullerides. If higher sub bands 
are included, one may argue that Ep is much larger, but 
the conclusion that Migdal's theorem is violated remains 
correct (Gunnarsson et al., 1994b). 

The lowest order vertex corrections have been studied 
by Grabowski and Sham (1984), who calculated these 
corrections approximately, finding that they can strongly 
reduce T c . More recently, the issue has been studied ex- 
tensively by Pietronero and coworkers (Pietronero and 
Strassler, 1992; Pietronero et al, 1995; Grimaldi et al, 
1995ab). They introduced a theory which is nonpertur- 
bative in A and perturbative in mX, where m = ujph/Ep. 



14 shows their results for T c and for the isotope effect 
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where M is the carbon mass. These results were calcu- 
lated assuming that there is a cut off in g-space at some 
critical Q c . It was argued (Grimaldi et al., 1995ab) that 
this cut off may simulate specific one-particle effects or 
correlation effects (Kulic and Zheyer, 1994). Depending 
on the value of Q c , the theory may predict a large en- 
hancement of T c or no enhancement at all. The value of 
the isotope effect can also vary substantially, and even be 
larger than 0.5. In the limit where Migdal's theorem is 
valid, this theory gives a = 0.5, since no Coulomb effects 
are included (fi* — 0). 
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FIG. 12. The zeroth order vertex and the first order vertex 
correction. (After Pietronero et al, 1995) 

First the lowest order vertex correction was calculated 
(see Fig. 12) for an electron gas like model. The result 
is shown in Fig. 13 for m = 1 and as a function of lo 
and Q = q/(2kp), where u and q are the Matsubara fre- 
quency and wave vector of the phonon line and kp is the 
Fermi wave vector. The energy of the incoming electron 
was set to zero, and w = uj ph is the energy of the Ein- 
stein phonon. This should be compared with the zeroth 
order contribution, which is unity. The vertex correc- 
tion has a very strong Q dependence. For small Q ^ 
and small u> it is large and negative, while it is positive 
elsewhere. This suggests that the influence on T c may 
depend strongly on the system. An Eliashberg-like the- 
ory was then developed, including the lowest order vertex 
correction and two "crossed" phonon lines in the scatter- 
ing between two electrons (Grimaldi et al.. 1995ab). Fig. 



.0.6' ■ ■ ■ 1 ' 

0.5 1 1.5 2 2.5 3 

CO / to 

o 

FIG. 13. Behavior of the vertex correction for different val- 
ues of Q = <j/(2/cf). (After Pietronero et al, 1995). 

As described in Sec. Vi. Takada (1993) also considered 
corrections to Migdal's theorem, and found that these en- 
hance T c for moderate values of the phonon frequencies. 
Ikeda et al. (1992) showed that Migdal's theorem be- 
comes valid again when the phonon energies are much 
larger than the electronic energies. Kostur and Mitrovic 

(1993) studied the vertex corrections with the empha- 
sis on two-dimensional systems and Krishnamurthy et al. 

(1994) and Cappelluti and Pietronero (1996) studied sys- 
tems with a van Hove singularity. Kostur and Mitrovic 
(1994) studied the corrections to T c from vertex correc- 
tions for several model interactions. For a BCS-type of 
interaction and an isotropic electron-phonon interaction, 
they found that vertex corrections reduce T c , while for a 
highly anisotropic interaction, inspired by the cuprate su- 
perconductors, they found that T c can also be enhanced 
by the vertex corrections. Asai and Kawaguchi (1992) ar- 
gued that nonadiabatic effect could enhance the coupling 
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by 10-30 %. 



B. Jahn-Teller effects beyond Migdal-Eliashberg treatment 



The violation of Migdal's theorem is also expected to 
play an important role for the photoemission spectrum. 
According to Migdal's theorem, the spectrum should just 
show single phonon satellites. Due to various broadening 
effects, this cannot be checked experimentally for A 3 Ceo- 
Theoretical calculations suggest, however, that multiple 
satellites develop (Gunnarsson et al, 1994ab). This in- 
volves a substantial transfer of spectral weight to higher 
binding energies, and the coupling to phonons and plas- 
mons may explain (Knupfcr et al., 1993) the large ob- 
served width of the t\ u band in photoemission (Chen et 
al., 1991). . The work above illustrates that the violation 
of Migdal's theorem could have important effects for the 
fullerides, but the magnitude, and perhaps even the sign, 
of the corrections seem to be not fully established. 
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FIG. 14. (a) T c as a function of rn = co p h/EF and the cut 
off momentum Q c for A = 0.5. (b) Isotope effect exponent a 
as a function of m and Q c . (After Grimaldi et al, 1995). 



The three-fold degeneracy of the ii„-level is lifted when 
a vibration of H 9 symmetry is excited. This makes Ceo 
molecules with a partly filled t\ u level unstable to a Jahn- 
Teller distortion (Jahn and Teller, 1937). It has been 
shown (Ihm, 1994; Auerbach et al., 1994; Manini et al, 
1994) that this Jahn-Teller system has a Berry phase 
(Berry, 1984), which imposes selection rules on the al- 
lowed transitions. These Jahn-Teller modes are included 
in the Migdal-Eliashberg theory under the assumption 
that the band width is much larger than the phonon 
frequency. Below we discuss how in the opposite limit, 
the molecular limit, some new interesting effects appear. 
These effects may possibly be relevant also for doped Ceo 
compounds, although the band width may be too large 
for these systems to allow the survival of the molecular 
properties. 

It is interesting to study the energy of Cg^~ as a func- 
tion of the electron-phonon coupling strength. A Jahn- 
Teller energy is defined as 
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where the sum is over the eight H 9 modes and g v and oj v 
are the coupling strength and energy, respectively, of the 
vth mode. In the strong-coupling limit, the energy can 
then be written as (Lannoo et al., 1991) 



E bc (n) =C{n)E JT , 



(31) 



where C(n) = 1 (n = 1,5), 4 (n = 2,4) or 3 (n = 3). 
It was emphasized by Yabana and Bertsch (1992), Auer- 
bach et al. (1994) and Manini et al. (1994) that in the 
weak-coupling limit the prcfactor in the energy expres- 
sion is a factor 5/2 larger 



E W °(n) = °-C{n)E JT , 



(32) 



due to the Jahn-Teller effect, while for a mode of A g 
symmetry the two expressions (31) and (32) are identi- 
cal. This result is valid for a free molecule or a solid 
with zero band width, while the Migdal-Eliashberg the- 
ory is valid in the opposite limit of a large band width. 
Auerbach et al. (1994) and Manini et al. (1994) ob- 
served that in the zero band width, weak-coupling limit, 
the attractive interaction between two electrons due the 
the electron-phonon interaction for C 6 q , 
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is a factor of three larger than in the Eliashberg the- 
ory. This could potentially have important implications 
for the superconductivity. We notice, however, that for 
realistic estimates of the coupling constants, the system 
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is not in the weak-coupling limit and may actually be 
closer to the strong-coupling limit (Gunnarsson, 1995). 
The prefactor in Eq. (32, 33) is then strongly reduced 
already in the zero band width limit. 

VII. CALCULATED PROPERTIES 
A. Transition temperature 

It was very early pointed out that the estimated values 
of the electron-phonon coupling A to the intramolecular 
modes are of the right order of magnitude to explain the 
transition temperature T c (Varma et al, 1991; Schluter 
et al, 1992; Mazin et al, 1992). Since the calculated 
values of A are only intermediate, it is important for the 
explanation of the rather high values of T c that the in- 
tramolecular phonon frequency are quite high. For values 
of A in the upper range of the estimates, the rather large 
values of T c can be understood even if it is assumed that 
the inefficient renormalization of /z* by retardation effects 
in the fullerides leads to substantially larger /i* than for 
conventional superconductors. For instance, Gunnars- 
son et al. (1995) solved the Eliashberg equation, using 
the coupling constants deduced from photoemission for 
a free Cg molecules (see Sec. III. A. 4), a t\ u band with 
the width 1/2 eV and the density of states at the Fermi 
energy deduced from NMR (see, Sec. III. A. 4). It was 
found that to reproduce the values of T c for K3C60 and 
RbaCgo, had to be assumed to be 0.6, i.e., even larger 
than the estimate in Sec. V.B. We note that the value 

0. 6 used in the Eliashberg equation corresponds to the 
value 0.4 in the McMillan formula, due to renormaliza- 
tion effects inside the t\ u band. From the discussion in 
the previous sections it is, however, clear that the un- 
certainties in the estimates of A and /i* are substantial, 
and the use of the Migdal-Eliashberg theory introduces 
further uncertainties. 

Zheng and Bennemann (1992) studied effects of the 
phonon energies being comparable to the band width. 
They used the Eliashberg theory, i.e., Migdal's theorem 
was implicitly assumed to be valid, but they did not as- 
sume the band width to be large in the solution of the 
Eliashberg theory. As a result they found that the quasi- 
particle weight was less renormalized than is normally 
assumed (the weight was larger than 1/(1 + A)). As a re- 
sult their calculated value of T c was larger than according 
to, e.g., the McMillan (1968) formula. 

1. Lattice parameter dependence 

As discussed in Sec. II. A, for the A3C60 superconduc- 
tors, except CS3C60, T c increases with the lattice parame- 
ters. The explanation of this behavior was provided very 
early (Fleming et al, 1991; Zhou et al, 1992; Varma 



et al, 1991; Schluter et al, 1992; Mazin et al, 1992). 
As discussed in Sec. III. A. 2, for the coupling to the in- 
tramolecular phonons A = VN(0) can to a very good 
approximation be factorized in an intramolecular part V 
and the density of states N(0) (Lannoo et al, 1991). The 
intramolecular part V can be expected to be essentially 
independent of the lattice parameter a, while N(0) grows 
with a. This follows, since as a increases, the overlap be- 
tween the molecules is reduced, leading to a reduced band 
width and an increased density of states. Fig. 15 shows 
a calculation of T c based on a tight-binding calculation 
of the N(0), assuming a power law reduction of the hop- 
ping matrix elements with the separation between the 
closest atoms on neighboring molecules (Schluter et al, 
1992b). The exponent of the power law was adjusted to 
experiment. 
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FIG. 15. Experimental (circles and crosses) and theoretical 
(full line) variation of T c with the shortest distance d between 
two closest carbon atoms on neighboring C60 molecules. The 
density of states was assumed to have a power law depen- 
dence on d, N(0) ~ d~ n , with n = 2.7 adjusted to experi- 
ment. (After Schluter et al, (1992b)) 

The very strong variation of T c with the lattice param- 
eter a in Na 2 AC6o (see Fig. 3) suggests a strong variation 
of N(0) with a in this system (see, e.g., Yildirim et al, 
1995). Indeed, Maniwa et al. (1995) deduced N(0) from 
the NMR relaxation rate for Na 2 AC 60 (A=K, Rb, Cs) 
at ambient pressure and found a much stronger depen- 
dence than for the Fm3m structure, although possibly 
not quite strong enough to fully explain the strong varia- 
tion of T c . Simple tight-binding calculations including 60 
"radial" p-orbital per Ceo molecule did, however, not give 
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a very different a-dependence for the density of states for 
the two structures (Gunnarsson, 1993). It has instead 
been suggested that Na2KC6o and Na2RbCgo may have 
structural distortions, and that these distortions possi- 
bly could reduce N(0) much more than the change in 
the lattice parameter (Maniwa, et al., 1995). Since the 
Fermi energy falls fairly close to a maximum in the den- 
sity of states in the undistorted Pa3 structure (Satpathy 
et al., 1992), such a reduction is possible. An intermedi- 
ate a-dependence was observed in pressure experiments 
for Na 2 CsC6o (Mizuki, 1994), possibly consistent with 
the observation that the structure may also distort un- 
der pressure (Zhu, 1995; Fischer, 1996). Whether N(0) 
in the ideal structure has a strong dependence on a or 
the strong dependence could be explained by other fac- 
tors remains an interesting issue. For instance, it has 
also been proposed that paramagnetic impurities may 
suppress T c in Na 2 KCeo and Na 2 RbCeo (Tanigaki and 
Prassides, 1995). 



2. Doping dependence 

As discussed in Sec. II. A, Yildirim et al. (1995, 
1996ab) have found a strong doping (n) dependence for 
T c , by studying Na2Cs x Ceo (n < 3) and Rba-^Ba^Cgo 
(n > 3). They found a maximum for T c close to half- 
filling (n = 3) and rapid decrease as n was reduced or 
increased. The reduction for n > 3 may be explained by 
considering the density of states for a system with mero- 
hedral disorder (Yildirim et al., 1995). The density of 
states shows a steady reduction with energy, except at 
the bottom of the band (Gelfand and Lu, 1992). As the 
filling n is increased, the density of states at the Fermi 
energy N(0) then drops, and A is reduced, which may 
explain the drop in T c . It appears, however, that there is 
no such explanation for n < 3. A tight-binding calcula- 
tion for the appropriate structure (Pa3) shows an increase 
in the density of states as the energy is lowered (Satpa- 
thy et al., 1992). The situation is changed slightly if 
the potential from the alkali atoms is taken into account 
(Gunnarsson, 1996). The nearest alkali neighbors of a 
given Ceo molecule set up a potential with I = 4, which 
cannot be efficiently screened by the Cgn molecule, and 
which modifies the density of states. If it is taken into ac- 
count, however, that the density of states is sampled over 
a finite energy, due to the rather large phonon energies, 
the basic tendency of the density of states to grow when 
the energy is reduced remains. One might have expected 
the potential of the vacancies to strongly influence the 
density of states. The potential from the vacancies has, 
however, mainly monopole or dipole character relative to 
a neighboring Ceo molecules, and it is therefore efficiently 
screened. Thus the drop in T c as n is reduced appears 
not to be understood. 



B. Isotope effect 

As discussed in Sec. II. C different experiment give 
widely different results for the isotope effect (T c ~ M~ a ), 
with a even being larger than unity. In the case of almost 
complete substitution a = 0.30 was found for K 3 Ceo and 
Rb 3 C 60 (Chen and Liebcr, 1992). 

It was pointed out early that a reasonable value of 
a <~ 0.3 can be obtained if the mechanism for super- 
conductivity is assumed to be the interaction with the 
intramolecular phonons (Schluter, 1992b). For instance, 
from the McMillan formula (1), the prefactor uj p h gives 
a = 0.5 and the renormalization of /i* as defined in Eq. 
(17) gives a reduction of a. Thus a is reduced to about 
0.3 if it is assumed that [i* <~ 0.2 and ui p h ~ 1000 cm -1 . 
We note, however, that the McMillan formula is not very 
accurate for calculating the isotope effect. The isotope 
effect was also calculated in the model discussed in the 
previous (VILA) section using the Eliashbcrg theory and 
the coupling constants derived from photoemission (Gun- 
narsson et al, 1995). The values of a were then found to 
be 0.32 and 0.37 for K 3 Ceo and RbsCen, respectively. 

Dcaven and Rokhsar (1993) addressed the large spread 
in the experimentally observed isotope effect a for differ- 
ent isotope substitions. They considered a model with 
harmonic phonons and a linear electron phonon coupling. 
Within the Migdal-Eliashbcrg framework, they showed 
that the anomolous isotope effect cannot be explained. 
In particular, they showed that A is independent of the 
isotope distribution. Anharmonic effects have been con- 
sidered by, e.g., Crespi et al. (1991) and Inada and Nasu 
(1992). These effects could in principle lead to a > 0.5. 

The isotope effect has traditionally been considered as 
a sign of an electron-phonon mechanism. It has, however, 
been pointed out that even for an electronic mechanism, 
there is an isotope effect due to the change of the elec- 
tronic structure by zero-point vibrations (Chakravarty et 
at, 1992a, Ashcroft and Cyrot, 1993). Ashcroft and Cy- 
rot (1993) emphasized the opposing effects of the transla- 
tional phonons and the librations on the electronic struc- 
ture. They argued that an increase in the carbon mass 
would tend to reduce the hopping integral due to the 
translational vibrations but increase it due to the libra- 
tions. They found that depending on the precise assump- 
tions about the parameters, a could take a broad range of 
values, including values larger than 0.5. In view of this it 
is not clear how much is proven by the reasonable values 
for a found from the intramolecular phonons alone. 



C. Reduced gap 

The superconductivity gap A is of interest, since it may 
give indications about strong-coupling effects. Gunnars- 
son et al, (1995) calculated the gap using the Eliashbcrg 



27 



equation in the model discussed above, using the cou- 
plings deduced from photoemission. They obtained the 
values 2A/T C = 3.59 and 3.66 for K 3 C 60 and Rb 3 C 60 , 
respectively. This is close to the BCS value 3.52 and in 
agreement with most recent measurements. 

The much larger values found in point contact tunnel- 
ing measurements have stimulated a substantial amount 
of work. Mazin et al. (1993b) studied a model with 
coupling to both intramolecular phonons (with frequen- 
cies around 1000 cm -1 ) and to low- lying modes with fre- 
quencies of the order 40 cm -1 . The latter modes could 
be librations, intermolecular or alkali modes. Due to the 
huge difference between the frequencies, this model shows 
very interesting properties. The low-lying modes hardly 
influence T c but have a strong influence on the reduced 
gap 2A/T C . Thus Mazin et al. (1993b) could use typical 
couplings to the intramolecular modes (Ai ~ 0.5) to pro- 
duce the measured values of T c , but at the same time use 
strong coupling to the low-lying modes (A2 ~ 2.7) to pro- 
duce a large reduced gap (2A/T C ~ 5). More recent work 
has suggested that the coupling to the low-lying modes 
is much weaker (see Sec. III.B, C and D) and that the 
reduced gap may be much smaller than found in point 
contact tunneling. The model of Mazin et al. (1993b) 
may therefore not be applicable to A3C60, but it is still 
a nice illustration of the unconventional properties that 
can be obtained if bosons of very different energies couple 
to the electrons. 

Mele and Erwin (1993) considered the detailed Fermi 
surface for a system with ordered molecules. The stud- 
ied an axial, nodeless order parameter with <i-symmctry. 
This led to the ratio 2A/T C ~ 4.6. The temperature de- 
pendence and the quasi-particle spectrum were also non- 
BCS. Since the solution, however, depends on specific 
features of the Fermi surface, it is not clear if it survives 
for the disordered system. 

Mele et al. (1994) studied the effects of orienta- 
tional disorder. They found that this introduces some 
gap anisotropy. Furthermore they obtained some reduc- 
tion of both T c and A in such a way that the reduced 
gap 2A/T C was increased to values of the order 3.8-4.0. 
They therefore concluded that disorder can lead to an 
increased gap ratio even within a weak-coupling theory. 

D. Other properties 

Erwin and Pickett (1991) calculated the Fermi veloc- 
ity vp for a system with ordered Ceo molecules, using 
the local density approximation. From up = 1.8 10 7 cm 
s _1 they deduced the Drude plasma energy Hil = 1.2 eV 
and the clean limit London penetration depth A = 1600 
A. Considering that A3C60 is in the dirty limit, they 
corrected this result by a factor (1 + ^/l) 1 ^ 2 , where £ is 
the coherence length and I is the mean free path. Us- 
ing £ = 26 A and I = 10 A, they obtained A - 3000 A, 



in reasonable agreement with experiment (the scattering 
between different experiments is, however, large). 



VIII. CONCLUDING REMARKS 

It is widely believed that the A3C60 (A=K, Rb) are 
s-wave BCS superconductors, driven by the coupling to 
the intramolecular H g phonons. The work reviewed here 
is consistent with and supports such a picture, although 
there is no conclusive evidence that this picture is correct 
or that an electronic mechanism is excluded. Accurate 
and detailed tunneling data would be very interesting in 
this context. 

It appears to be understood why A 3 Cgo are metals 
and not Mott-Hubbard insulators, in spite of the large 
Coulomb interaction. Experimental and theoretical es- 
timates of the electron-phonon coupling A give values of 
about 0.5-1. The spread in these estimates indicate that 
the coupling strength is still not very accurately known. 
For the other important parameter determining T c , the 
Coulomb pseudopotcntial there are arguments that 
/1* is not strongly renormalized by retardation effects due 
to scattering of the electrons into sub bands other than 
the t\ u band. This, however, leaves interesting questions 
about the screening and retardation effects on ji* within 
the ti u band. Theoretical calculations of T c based on es- 
timates of A and fi* are consistent with experiment and 
support an electron-phonon picture. 

The lattice parameter dependence of T c seems to be 
understood for A3C60 in the Fm3m structure but not for 
Na 2 AC6o in the Pa3 structure. The doping dependence is 
also not well understood. The BCS like superconducting 
gap seen in most experiment, follows from the picture of 
the superconductivity being driven by the H g intramolec- 
ular phonons, while the isotope effect still gives rise to 
some interesting experimental and theoretical questions. 

We have emphasized that because several energy scales 
are comparable, the theoretical treatment of the ful- 
lerides is particularly demanding but also unusually in- 
teresting. Thus one may expect Migdal's theorem to be 
violated and many-body effects due to the Coulomb inter- 
action to be important. The electron-phonon interaction 
should also be important for many properties. On the 
other hand, the molecular solid character leads to cer- 
tain simplifications, for instance the electronic properties 
are expected to be dominated by the three-fold degen- 
erate ti u orbital. Due to the unusual parameter range, 
many interesting general issues are raised in a particular 
clear-cut way for the doped fullerides, and their solutions 
should have implications in many other fields. 
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